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Abstract 

We review a class of non-topological defects in the standard electroweak model, and their 
implications. Starting with the semilocal string, which provides a counterexample to many 
well known properties of topological vortices, we discuss electroweak strings and their sta- 
bility with and without external influences such as magnetic fields. Other known properties 
of electroweak strings and monopoles are described in some detail and their potential rele- 
vance to future particle accelerator experiments and to baryon number violating processes is 
considered. We also review recent progress on the cosmology of electroweak defects and the 
connection with superfluid helium, where some of the effects discussed here could possibly 
be tested. 
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1 Introduction 



In a classic paper from 1977 ||9i|, a decade after the SU{2)l x U{1)y model of electroweak 
interactions had been proposed |}49(| , Nambu made the observation that, while the Glashow- 
Salam- Weinberg (GSW) model does not admit isolated, regular magnetic monopoles, there 
could be monopole-antimonopole pairs joined by short segments of a vortex carrying Z- 
magnetic field (a Z- string). The monopole and antimonopole would tend to annihilate but, 
he argued, longitudinal collapse could be stopped by rotation. He dubbed these configu- 
rations dum&eH^ and estimated their mass at a few TeV. A number of papers advocating 
other, related, soliton-type solutions Q in the same energy range followed |Q, but the lack 
of topological stability led to the idea finally being abandoned during the eighties. 



Several years later, and completely independently, it was observed that the coexistence of 
global and gauge symmetries can lead to stable non-topological strings called "semilocal 
strings" |121| in the sva^ 9^ = 1 limit of the GSW model that Nambu had considered. 
Shortly afterwards it was proved that Z-strings were stable near this limit |117|, and the 
whole subject made a comeback. This report is a review of the current status of research 
on electroweak strings. 



Apart from the possibility that electroweak strings may be the first solitons to be observed in 
the standard model, there are two interesting consequences of the study of electroweak and 
semilocal strings. One is the unexpected connection with baryon number and sphalerons. 
The other is a deeper understanding of the connection between the topology of the vacuum 
manifold (the set of ground states of a classical field theory) and the existence of stable 
non-dissipative configurations, in particular when global and local symmetries are involved 
simultaneously. 



In these pages we assume a level of familiarity with the general theory and basic properties of 
topological defects, in particular with the homotopy classification. There are some excellent 
reviews on this subject in the literature to which we refer the reader |l05[ |2|, On the 
other hand, electroweak and semilocal strings are non-topological defects, and this forces 
us to take a slightly different point of view from most of the existing literature. Emphasis 
on stability properties is mandatory, since one cannot be sure from the start whether these 
defects will actually form. With very few exceptions, this requires an analysis on a case by 
case basis. 



Following the discussion in |Q, one should begin with the definition of dissipative config- 
urations. Consider a classical field theory with energy density Too > such that Too = 
everywhere for the ground states (or "vacua" ) of the theory. A solution of a classical field 
theory is said to be dissipative if 

lim maxxToo(x, i) = (1) 



We will consider theories with spontaneous symmetry breaking from a Lie group G (which 
we assume to be finite-dimensional and compact) to a subgroup H] the space V of ground 
states of the theory is usually called the vacuum manifold and, in the absence of accidental 
degeneracy, is given by V = G/H . 

^or monopolia, after analogous configurations in superfluid helium 

^One example, outside the scope of the present review, are so-called vorticons, 

proposed by Huang and Tipton, which are closed loops of string with one quantum 

of Z boson trapped inside. 
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The classification of topological defects is based on the homotopy properties of the vacuum 
manifold. If the vacuum manifold contains non-contractible n-spheres then field configura- 
tions in n + 1 spatial dimensions whose asymptotic values as r — > cxd "wrap around" those 
spheres are necessarily non-dissipative, since continuity of the scalar field guarantees that, 
at all times, at least in one point in space the scalar potential (and thus the energy) will be 
non-zero. The region in space where energy is localized is referred to as a topological defect. 
Field configurations whose asymptotic values are in the same homotopy class are said to be 
in the same topological sector or to have the same winding number. 

In three spatial dimensions, it is customary to use the names monopole, string and domain 
wall to refer to defects that are pointlike, one-dimensional or two-dimensional respectively. 
Thus, one can have topological domain walls only if 7ro(V) ^ 1, topological strings only 
if 7ri(V) ^ 1 and topological monopoles only if 7T2{V) ^ 1. Besides, defects in different 
topological sectors cannot be deformed into each other without introducing singularities or 
supplying an infinite amount of energy. This is the origin of the homotopy classification 
of topological defects. We should point out that the topological classification of textures 
based on 7r3(V) has a very different character, and will not concern us here; in particular, 
configurations from different topological sectors can be continuously deformed into each 
other with a finite cost in energy. In general, textures unwind until they reach the vacuum 
sector and therefore they are dissipative. 

It is well known, although not always sufficiently stressed, that the precise relationship 
between the topology of the vacuum and the existence of stable defects is subtle. First of 
all, note that a trivial topology of the vacuum manifold does not imply the non-existence 
of stable defects. 



Secondly, we have said that a non-trivial homotopy of the vacuum manifold can result in 
non-dissipative solutions but, in general, these solutions need not be time independent nor 
stable to small perturbations. One exception is the field theory of a single scalar field in 1+1 
dimensions, where a disconnected vacuum manifold (i.e. one with 7ro(V) ^ 1) is sufficient 
to prove the existence of time independent, classically stable "kink" solutions |3^. But 
this is not the norm. The 0(3) model, for instance, has topological global monopoles [|l5j 
which are time independent, but they are unstable to angular collapse even in the lowest 
non-trivial winding sector |^o| . 

It turns out that the situation is particularly subtle in theories where there are global and 
gauge symmetries involved simultaneously. The prototype example is the semilocal string, 
described in section ^j. In the semilocal string model, the classical dynamics is governed by 
a single parameter /3 = m^/m^ that measures the square of the ratio of the scalar mass, 
TOs, to the vector mass, ruy (this is the same parameter that distinguishes type I and type 
II superconductors). It turns out that: 

- when /3 > 1 the semilocal model provides a counterexample to the widespread belief that 
quantization of magnetic flux is tantamount to its localization, i.e., confinement. The vector 
boson is massive and we expect this to result in confinement of magnetic flux to regions 
of width given by the inverse vector mass. However, this is not the case! As pointed out 



by Hindmarsh te5| and Preskill |102 , this is a system where magnetic flux is topologically 



conserved and quantized, and there is a finite energy gap between the non-zero flux sectors 



^The names cosmic string and vortex are also common. Usually, "vortex" refers 
to the configuration in two spatial dimensions, and "string" to the corresponding 
configuration in three spatial dimensions; the adjective "cosmic" helps to distinguish 
them from the so-called fundamental strings or superstrings. 
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and the vacuum, and yet there are no stable vortices. 



- when /3 < 1 strings are stableQ even though the vacuum manifold is simply connected, 
'^lO^) — 1- Semilocal vortices with f3 < 1 are a remarkable example of a non-topological 
defect which is stable both perturbatively and to semiclassical tunnelling into the vacuum 
p3|- 



As a result, when the global symmetries of a semilocal model are gauged, dynamically 
stable non-topological solutions can still exist for certain ranges of parameters very close to 
stable semilocal limits. In the case of the standard electroweak model, for instance, strings 
are (classically) stable only when sin^ 6*^, « 1 and the mass of the Higgs is smaller than the 
mass of the Z boson. 



We begin with a description of the Glashow-Salam- Weinberg model, in order to set our 
notation and conventions, and a brief discussion of topological vortices (cosmic strings). 
It will be sufficient for our purposes to review cosmic strings in the Abelian Higgs model, 
with a special emphasis on those aspects that will be relevant to electroweak and semilocal 
strings. We should point out that these vortices were first considered in condensed matter 
by Abrikosov in the non-relativistic case, in connection with type II superconductors. 
Nielsen and Olesen were the first to consider them in the context of relativistic field theory, 
so we will follow a standard convention in high energy physics and refer to them as Nielsen- 
Olesen strings Q . 

Sections || to |^ are dedicated to semilocal and electroweak strings, and other embedded 
defects in the standard GSW model. Electroweak strings in extensions of the GSW model 
are discussed in section |[ 

In section 1^ the stability of straight, infinitely long electroweak strings is analysed in detail 
(in the absence of fermions). Sections ^ to ^ investigate fermionic superconductivity on 
the string, the effect of fermions on the string stability, and the scattering of fermions off 
electroweak strings. The surprising connection between strings and baryon number, and 
their relation to sphalerons, is described in sections || and Here we also discuss the 
possibility of string formation in particle accelerators (in the form of dumbells, as was 
suggested by Nambu in the seventies) and in the early universe. 

Finally, section ^ describes a condensed matter analog of electroweak strings in superfluid 
helium which may be used to test our ideas on vortex formation, fermion scattering and 
baryogenesis. 



A few comments are in order: 



• Unless otherwise stated we take spacetime to be flat, 3+1 dimensional Minkowski space; 
the gravitational properties of embedded strings are expected to be the same as those of 
Nielsen- Olesen strings and will not be considered here. A limited discussion of possible 
cosmological implications can be found in sections 3.5 and 9.4. 



• We concentrate on regular defects in the standard model of electroweak interactions. 
Certain extensions of the Glashow-Salam- Weinberg model are briefly considered in section 
6 but otherwise they are outside the scope of this review; the same is true of singular 
solutions. In particular, we do not discuss isolated monopoles in the GSW model p8[ |, 

^We want to stress that, contrary to what is often stated in the literature, the 
semilocal string with /? < 1 is absolutely stable, and not just metastable. 
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which are necessarily singular. 

• No family mixing effects are discussed in this review and we also ignore SU{5)c colour 
interactions, even though their physical affects are expected to be very interesting, in par- 
ticular in connection with baryon production by strings (see section ^ . 

• Our conventions are the following: spacetime has signature (+, — , — , — ). Planck's constant 
and the speed of light are set to one, h ~ c ^ 1. The notation (x) is shorthand for all 
spacetime coordinates (a;'^,a;*), i = 1,2,3; whenever the ^-coordinate is meant, it will be 
stated explicitly. We also use the notation {t,x). 

• Complex conjugation and hermitian conjugation are both indicated with the same symbol, 
(^), but it should be clear from the context which one is meant. For fermions, ip = V'W^j 
usual. Transposition is indicated with the symbol (^). 

• One final word of caution: a gauge field is a Lie Algebra valued one- form A — A^^dx^^ — 
A'l^^T'^dx^, but it is also customary to write it as a vector. In cylindrical coordinates 
(t, p, ip, z), A = Atdt + Apdp+ A^d(p + Azdz is often written A = Ati+App+{A^/p)(f + AzZ, 
In spherical coordinates, (t,r,d,(p), A — Atdt + A^dr + Agd9 + A^dtp is also written 
A — Aft + Arf + {Ag/r)6 + {A^/r sin 9)(p. We use both notations throughout. 



1.1 The Glashow-Salam- Weinberg model. 



In this section we set out our conventions, which mostly follow those of |27| . 

The standard (GSW) model of electroweak interactions is described by the Lagrangian 

L = Lb + + (2) 

families 

The first term describes the bosonic sector, comprising a neutral scalar 4>^ , a charged scalar 
0+, a massless photon A^, and three massive vector bosons, two of them charged [W^) 
and the neutral Z^. 

The last two terms describe the dynamics of the fermionic sector, which consists of the three 
families of quarks and leptons 













T 






t 









(3) 



1.1.1 The bosonic sector 

The bosonic sector describes an SU{2)l x U{\)y invariant theory with a scalar field $ in 
the fundamental representation of SU{2)l. It is described by the Lagrangian: 

Lb = Lw + Ly + -V{^) (4) 

with 

= - V.^'^W^^"" , a = 1,2, 3 

f (5) 



8 



where W^J" = d^,W^ - d^W^ + ge'^^'^W^^W;^ and Y^^ = a^F^ - 9^Yf, are the field 
strengths for the SU{2)l and U{1)y gauge fields respectively. Summation over repeated 
SU{2)l indices is understood, and there is no need to distinguish between upper and lower 



ones, e 



123 _ 



1. Also, 



\D,<f\^ 



(6) 
(7) 



where r" are the Pauli matrices, 



ri = 



1 

1 



r2 = 



-i 

1 



1 

- 1 



(8) 



from which one constructs the weak isospin generators = ^t'' satisfying [T°,r''] = 



The classical field equations of motion for the bosonic sector of the standard model of the 
electroweak interactions are (ignoring fermions): 

+ 2A _ $ ^ (9) 

D^W"'' = j'^ = ^-g [^V'^D^* - (1?'^$)^°$] (10) 

where D^W>''"' = ^^W^''"' + ge'^^^W^W'''"' . 

When the Higgs field <i> acquires a vacuum expectation value (VEV), the symmetry breaks 
from SU{2) l x U{1)y to U{l)em- In particle physics it is standard practice to work in unitary 
gauge and take the VEV of the Higgs to be ($-^) = 77(0, 1)/V2. In that case the unbroken 
U{1) subgroup, which describes electromagnetism, is generated by the charge operator 

and the two components of the Higgs doublet are charge eigenstates 




Y is the hypercharge operator, which acts on the Higgs like the 2x2 identity matrix. 
Its eigenvalue on the various matter fields can be read off from the covariant derivatives 
Df, = df,- igWl^T" - ig'Y^,{Y/2) which are listed explicitly in equations (|) and (||)-(||). 

In unitary gauge, the Z and A fields are defined as 

= cos e.^Wl - sin e.u,Y^, , = sin 61^1^^ + cos O^Y^, , (14) 

and = (W^ =F iW^)/\/2 are the W bosons. The weak mixing angle 9w is given by 
tan^^ = electric charge is e = sin^^ cos^m with = {g"^ + g'^Y^^ ■ 

However, unitary gauge is not the most convenient choice in the presence of topological 
defects, where it is often singular. Here we shall need a more general definition in terms of 
an arbitrary Higgs configuration $(a;): 

= cos n''{x)W^ - sin 9^ , A^ = sin 9.^ n"(a;)VK; + cos 6^ Y^ , (15) 
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where 

- $t(x)$(x) ■ ^ ^ 

is a unit vector by virtue of the Fierz identity X]a(*^^'''"*^)'^ = (<i>'f$)^ . In what follows we 
omit writing the x-dependence of explicitly. Note that n° is ill defined when <i> = 0, so 
in particular at the defect cores. 

The generators associated with the photon and the Z-boson are, respectively, 

Q = n^T'' + Y/2 , Tz = cos^ Q^n^T" - sin^ oj^^ n'^T" ~ sit? 9^Q , (17) 

while the generators associated with the (charged) W bosons are determined, up to a phase, 
by the conditions 

[Q, r±] = ±r± [T+, T-] = n^T" = Tz + sir? 6^Q , (T+)'f = T" (18) 

(note that, if = (0, 0, 1) as is the case in unitary gauge, one would take = {T^ ± 

There are several different choices for defining the electromagnetic field strength but, fol- 
lowing Nambu, we choose: 

= sine^ n^WI^^ + cosByj F^^ (19) 

where, and Y^^, are field strengths. The different choices for the definition of the field 
strength agree in the region where 1?^$ — where is the covariant derivative operator; 
in particular this is different from the well known 't Hooft definition which is standard for 



monopoles [|6l|. (For a recent discussion of the various choices see, e.g. ^ IISD- And 
the combination of SU{2) and U{1) field strengths orthogonal to is defined to be the 
Z field strength: 

Z^, = COS0U, n'^Wl, - sin9^ Y^, . (20) 



1.1.2 The fermionic sector 



The fermionic Lagrangian is given by a sum over families plus family mixing terms {Lf„i). 
Family mixing effects are outside the scope of this review, and we will not consider them 
any further. Each family includes lepton and quark sectors 



Lf — Li + Lq 

which for, say, the first family are 



where 



- Gd {u,d)L 

- Gu {u, d)L 



t'R + any 



(21) 



(22) 



(23) 



where 0* and (j) are the complex conjugates of (jp and 0+ respectively, /i, Gd and Gu are 
Yukawa couplings. The indices L and R refer to left- and right-handed components and. 
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rather than hst their charges under the various transformations, we give here all covariant 
derivatives explicitly: 



D^dR = (d^ + '-^Y^)dR 



(24) 

(25) 
(26) 

(27) 
(28) 



• One final comment: 



Electroweak strings are non-topological and their stability turns out to depend on the values 
of the parameters in the model. In this paper we will consider the electric charge e, Yukawa 
couplings and the VEV of the Higgs, ri/\/2, to be given by their measured values, but the 
results of the stability analysis will be given as a function of the parameters sin^ 6^ and 
P = {mji /mzY' (the ratio of the Higgs mass to the Z mass squared); we remind the reader 
that sin^ 6^ « 0.23, mz = gzVf^ = 91-2 GeV, mw = 5'7/2 = 80.41 GeV and current bounds 
on the Higgs mass rriH = V^Xrj are rriH > 77.5 GeV. 



2 Reviev^r of Nielsen-Olesen topological strings 



We begin by reviewing Nielsen-Olesen (NO) vortices in the Abelian Higgs model, with 
emphasis on those aspects that are relevant to the study of electroweak strings. More detailed 
information can be found in existing reviews |105|. 



2.1 The Abelian Higgs model 



The theory contains a complex scalar field <i> and a J7(l) gauge field which becomes massive 
through the Higgs mechanism. By analogy with the GSW model, we will call this field Y^. 
The action is 



S 



d^x 



1 V Vt^'^ 
4 -[ /j;/ ^ 



(29) 



where = d^ — iqY^ is the C/(l)-covariant derivative, and Y^^ = d^Y^ — d^Y^^ is the U{\) 
field strength. The theory is invariant under U{\) gauge transformations: 

$(:;,) ^ e*«^(")$(a;) = ^x) , Y^{x) Y^{x) + d^x{x) = Y^,{x) (30) 
which give D^$(a;) -> D^^{x) = e*«x("')D^<I>. 



11 



The equations of motion derived from this Lagrangian are: 



Df,D^'^ + 2A(|$|2 - = 



(31) 



Before we proceed any further, we should point out that, up to an overall scale, the classical 
dynamics of the Abelian Higgs model is governed by a single parameter, (3 = 2X/q^, the 
(square of the) ratio of the scalar mass to the vector mass ^. The action (|9|) contains three 
parameters, (A, 77,(7), which combine into the scalar mass \f2\r\ = rris = /j^, the vector 



mass qr\ 



the Higgs, ri/^/2. The rescaling 



$(x) = -^$(x) 



^, and an overall energy scale given by the vacuum expectation value of 



V2 



— X 
qrj 



Y - —Y 



changes the action to 



(32) 



(33) 



where now = d^ — lY^^ and we have omitted hats throughout for simplicity. In physical 
terms this corresponds to taking Z„ as the unit of length (up to a factor of ^/2) and absorbing 
the C/(l) charge q into the definition of the gauge field, thus 



< m > 



E 



E 



< l$p > 



(34) 



The energy associated with (03) is 



(35) 



where the electric and magnetic fields are given by Foi = Ei and Fij — €ijkB^ respectively 
{i,j,k = 1,2,3). Modulo gauge transformations, the ground states are given by = 0, 
$ = r]e^'~' /^/2, where C is constant. Thus, the vacuum manifold is the circle 



V = {$ e (C I - y = 0} 51 



(36) 



A necessary condition for a configuration to have finite energy is that the asymptotic scalar 
field configuration as r ^ 00 must lie entirely in the vacuum manifold. Also, -D^$ must tend 
to zero, and this condition means that scalar fields at neighbouring points must be related 
by an infinitesimal gauge transformation. Finally, the gauge field strengths must also vanish 
asymptotically. Note that, in the Abelian Higgs model, the last condition follows from the 
second, since — [D^,D,y]^ = —iqYfj,,^^ implies Y^jy — But this need not be the case 
when the Abelian Higgs model is embedded in a larger model. 

Vanishing of the covariant derivative term implies that, at large r, the asymptotic configu- 
ration ^{x) must lie on a gauge orbit; 

^(x) = 5(x)$o , where g{x) G G and $0 G V . (37) 

^ 13 is also the parameter that distinguishes superconductors or type I {(3 < 1) 
from type II (/? > 1) 
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where $o is a reference point in V. Note that, since all symmetries are gauge symmetries, 
the set of points that can be reached from <I>o through a gauge transformation (the gauge 
orbit of $o) spans the entire vacuum manifold. Thus, V = G/H — Giocai/^^iocai, where 
Giocai indicates the group of gauge - i.e. local - symmetries. On the other hand the spaces 
V and GiocaiZ-ff local need not coincide in models with both local and global symmetries, and 
this fact will be particularly relevant in the discussion of semilocal strings. 



2.2 Nielsen- Olesen vortices 



In what follows we use cylindrical coordinates {t,p,ip,z). We are interested in a static, 
cylindrically symmetric configuration corresponding to an infinite, straight string along the 
z-axis. 

The ansatz of Nielsen and Olesen for a string with winding number n is 

$ = -^/(p)e^"'^ , qY^^nvip) , Yp^Yt^Y.^O (38) 

(that is, Y = v{p) dip or Y ~ tp v{p)/ p ), with boundary conditions 



/(O) = ^;(0) = , /(p) ^ 1 , v{p) ^1 as p ^ oo. (39) 

Note that, since Y^ — Yt — 0, and all other fields are independent of t and z, the electric 
field is zero, and the only surviving component of the magnetic field B is in the z direction. 

Substituting this ansatz into the equations of motion we obtain the equations that the 
functions / and v must satisfy: 



/"(p) + IM^ _ + xvHi f{pr)f{p) - 



In what follows we will denote the solutions to the system ( ^ , ^9| ) by /jvo a-nd wato; they 
are not known analytically, but have been determined numerically; for n — 1, (3 = 0.5, they 
have the profile in Fig. 0. 

At small p, the functions / and v behave as p" and p^ respectively; as p — > oo, they approach 
their asymptotic values exponentially with a width given by the inverse scalar mass, m^, 
and the inverse vector mass, rriy, respectively, if /3 < 4. For /3 > 4 the fall-off of both the 
scalar and the vector is controlled by the vector mass 

One case in which it is possible to find analytic expressions for the functions f^o and vno 
is in the limit n — > cxj Inside the core of a large n vortex, the functions / and v are 

/(p) = (±.m,my) * e-'^"^"-'''/^ , v{p) = ^m^m^p^ (41) 

to leading order in 1/n, and the transition to their vacuum values is controlled by a first 
integral 5'(/, /', v, v') — const . Large n vortices behave like a conglomerate of "solid" n — 1 
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Figure 1: The functions /ato, vno for a string with winding number n = 1 (top panel) and 
n = 50 (bottom panel), for f3 = 2X/q^ = 0.5. The radial coordinate has been rescaled as in 
eq. (HI), p = qVP/V2. 
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vortices. The area scales as n, so the radius goes hke ^/nLo, where Lq = 2{^/rrr^rn^)^^ . The 
transition region between the core and asymptotic values of the fields is of the same width 
as for n = 1 vortices Fig. |l| shows the functions Jno, vnq for n = 50, (3 = 0.5 (note that 
for P > 1 these multiply winding solutions are unstable to separation into n = 1 vortices 
which repel one another). 



• Energy considerations: 

The energy per unit length of such configurations (static and z-independent) is therefore 



E= d 



j2. 



|A„$|' + ii3' + A($t<i.-!^)2 
where m,n — 1,2 and B = dmYn — dnY^ is the z-component of the magnetic field. 



(42) 



In order to have solutions with finite energy per unit length we must demand that, as 
p oo, -D^$, |$P — »7^/2 and y„i„ all go to zero faster than 1/p. 



The vacuum manifold (|3^) is a circle and strings form when the asymptotic field configu- 
ration of the scalar field winds around this circle. The important point here is that there 
is no way to extend a winding configuration inwards from p = oo to the entire xy plane 
continuously while remaining in the vacuum manifold. Continuity of the scalar field implies 
that it must have a zero somewhere in the xy plane. This happens even if the xy plane is 
deformed, and at all times, and in three dimensions one finds a continuous line of zeroes 
which signal the position of the string (a sheet in spacetime) . Note that the string can have 
no ends; it is cither infinitely long or a closed loop. 

The zeroes of the scalar field are forced by the non-zero topological degree of the map 

^ V 

(43) 

ip $(p = CX),(^) , 

usually called the winding number of the vortex; the resulting vortices are called topo- 
logical because they are labelled by non-trivial elements of the first homotopy group of 
the vacuum manifold (where non-trivial means "other than the identity element"). Thus, 
7ri(V) = 7ri(S'^) 1, is a necessary condition for the existence of topological vortices. Vor- 
tices whose asymptotic scalar field configurations are associated with the identity element 
of 7ri(V) are called non-topological. In particular, if V is simply connected, i.e. 7ri(V) = 1, 
one can only have non-topological vortices. 

A few comments are needed at this point: 



• Quantization of magnetic fiux: 

Recall that B is the ^-component of the magnetic field. The magnetic flux Fy through the 
ccy-plane is therefore 

Fy= f <fxB = [ Yoo-dl^ [ d^xd^ = — (44) 

and is quantized in units of 2'K/q. This is due to the fact that ^{p — oo, (p) — rje^i^^'^^ / \/2, 
D^^ — r\l \/2\iqd^x — i^Y^p] = and $ must be singlevalued, thus (7[x(27r) — x(0)] — ^ttti. 
The integer n is, again, the winding number of the vortex. 
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• Magnetic pressure: 

In an Abelian theory, the condition W ■ B = imphes that parallel magnetic field lines 
repel. A two-dimensional scale transformation x — > Xx where the magnetic field is reduced 
accordingly to keep the magnetic flux constant, B\ = hr'^B{x / K), reduces the magnetic 
energy J (Px _B^/2 by A^. What this means is that a tube of magnetic lines of area Sq 
can lower its energy by a factor of A'^ by spreading over an area A^5o. 

Note that later we will consider non- Abelian gauge symmetries, for which V • i? 7^ and the 
energy can also be lowered in a different way. In this case, one can think of the gauge fields 
as carrying a magnetic moment which couples to the "magnetic" field and, in the presence of 
a sufficiently intense magnetic field, the energy can be lowered by the spontaneous creation 
of gauge bosons. In the context of the electroweak model, this process is known as W- 
condensation and its relevance for electroweak strings is explained in section 0. 



• Meissner effect and symmetry restoration: 



In the Abelian Higgs model, as in a superconductor, it is energetically costly for magnetic 
fields to coexist with scalar fields in the broken symmetry phase. Superconductors exhibit 
the Meissner effect (the expulsion of external magnetic fields), but as the sample gets larger 
or the magnetic field more intense, symmetry restoration becomes energetically favourable. 
An example is the generation of Abrikosov lattices of vortices in type II superconductors, 
when the external magnetic field reaches a critical value. 



The same phenomenon occurs in the Abelian Higgs model. In a region where there is a 
concentration of magnetic fiux, the coupling term q^A^<^^ in the energy will tend to force 
the value of the scalar field towards zero (its value in the symmetric phase). This will 
be important to understand the formation of semilocal (and possibly electroweak) strings, 
where there is no topological protection for the vortices, during a phase transition (see 
section 3.5). The back reaction of the gauge fields on the scalars depends on the strength of 
the coupling constant q. When q is large (in a manner that will be made precise in section 
3.5) semilocal strings tend to form regardless of the topology of the vacuum manifold. 



2.3 Stability of Nielsen- Olesen vortices 



Given a solution to the classical equations of motion, there are typically two approaches 
to the question of stability. One is to consider the stability with respect to infinitesimal 
perturbations of the solution. If one can establish that no perturbation can lower the energy, 
then the solution is called classically stable. Small perturbations that do not alter the energy 
are called zero modes, and signal the existence of a family of configurations with the same 
energy as the solution whose stability we are investigating (e.g. because of an underlying 
symmetry). If one can guess an instability mode, this approach is very efficient in showing 
that a solution is unstable (by finding the instability mode explicitly) but it is usually much 
more cumbersome to prove stability; mathematically the problem reduces to an eigenvalue 
problem and one often has to resort to numerical methods. A stability analysis of this type 
for Nielscn-Olesen vortices has only been carried out recently by Goodband and Hindmarsh 
[ p2[ . An analysis of the stability of semilocal and electroweak strings can be found in later 
sections. 

A second approach, due to Bogomolnyi, consists in finding a lower bound for the energy 
in each topological sector and proving that the solution under consideration saturates this 
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bound. This immediately implies that the solution is stable, although it does not preclude 
the existence of zero modes or even of other configurations with the same energy to which 
the solution could tunnel semiclassically. We will now turn to Bogomolnyi's method in the 
case of Nielsen-Olesen vortices. 



• Bogomolnyi limit and bounds 
Consider the scalar gradients: 

{Di^)^Di'i> + {D2(S>)'' 02^^ = [(£>! + iD2)^]\Di + iDj)* - i[(Di$)tD2$ - (02^)'' Di^] 

(45) 

Note that the second term in the RHS of ( ^5|) is the curl of the current Ji — —i^^Di^, and 
that § J ■ dl tends to zero as p ^ oo for configurations with finite energy per unit length 
(because Di<^ must vanish faster than 1/p). Now use the identity [Di,£'2]$ — —iqFi2^ = 
—iqB^. to rewrite the energy per unit length as follows: 



E = I d^x 



d^3 



{Dl ± iD2)^^ + ^B^ ± qB<^>^^ + A _ 



\{Di±iD2)<i> 



B±q\ $t$ 



2 

{\-\q') \ 



±qL. j Bd^x 

(46) 

The last integral is the total magnetic flux, and we saw earlier that it has to be an integral 
multiple of 2Tr/q, so we can write, introducing /3 — 2X/q^, 



E = 2ni±n)^+J 



\{Di±iD2)'S? 



+ iq^(/3-l)|<f>T$-^ 



where the plus or minus signs are chosen so that the first term is positive, depending on 
the sign of the magnetic flux. 



Note that, if > 1 the energy is bounded below by 

E > ($t$)qFy 

where Fy is the magnetic flux. ^ 



(48) 



li f3 — 1, there are configurations that saturate this bound: those that satisfy the first order 
Bogomolnyi equations 

2 

(L>i±iD2)$ = , B±g($^$ - y) . (49) 



^ When /3 = 1, the masses of the scalar and the vector are equal, and the Abelian 
Higgs model can be made super symmetric. In general, bounds of the form (Energy) 
> (constant) x (flux) are called Bogomolnyi bounds, and their origin can be traced 
back to supersymmetry. 
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or, in terms of f{p) and v{p), 

/'(p) + (±n)^^^^/(p) = , (^±n)v'{p) + ^p{f{p)-l) = (50) 

p 2 



However, when /3 > 1 there does not exist a static solution with E = Tr\n\r]'^ since requiring, 
e.g., B + g($t$ - ?7^/2) = and {^^^ - ?7^/2) = simultaneously would imply B = 0, 
which is inconsistent with the condition on the total magnetic flux, J BcPx — 2Trn/q . This 
has an effect on the stability of higher winding vortices when /? > 1: if > 1 the solution 
breaks into n vortices each with a unit of magnetic flux | |2^ , which repel one another. 

If n = 1 there are stable static solutions, but with an energy higher than the Bogomolnyi 
bound. This is because the topology of the vacuum manifold forces a zero of the Higgs 
field, and then competition between magnetic and potential energy fixes the radius of the 
solution. The same argument shows that n ~ 1 strings are stable for every value of (3. One 
still has to worry about angular instabilities, but a careful analysis by shows there are 
none. 



The dynamics of multivortex solutions is governed by the fact that when f3 < 1 vortices at- 
tract, but with P > 1 they repel. This can be understood heuristically from the competition 
between magnetic pressure and the desire to minimise potential energy by having symmetry 
restoration in as small an area as possible. The width of the scalar vortex depends on the 
inverse mass of the Higgs, Ig, that of the magnetic flux tube depends on the inverse vector 
boson mass, 1^. If /3 < 1, have m„ > nis so 1^ < Ig (the radii of the scalar and vector 
tubes). The scalar tubes see each other first - they attract. Whereas if /3 > 1, the vector 
tubes see each other first - they repel. For /3 — 1 there is no net force between vortices, 
and there are static multivortex solutions for any n. In the Abelian Higgs case they were 
explicitly constructed by Taubes |Q and their scattering at low kinetic energies has been 
investigated using the geodesic approximation of Manton by Ruback [108| and, more 
recently, Samols [110|. For /3 < 1 Goodband and Hindmarsh |^ have found bound states 
of two n — I vortices oscillating about their centre of mass. 



3 Semilocal strings 



The semilocal model is obtained when we replace the complex scalar field in the Abelian 
Higgs model by an TV-component multiplet, while keeping only the overall phase gauged. In 
this section we will concentrate on = 2 because of its relationship to electroweak strings, 
but the generalisation to higher N is straightforward, and is discussed below. 



3.1 The model 



Consider a direct generalization of the Abelian Higgs model where the complex scalar field 

(01,02)- The action is 



is replaced by an SU{2) doublet $^ = ( 



s 



1 



(51) 
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where is the U{1) gauge potential and Y^i, = d^jYi, — d^Y^ its field strength. Note that 
this is just the scalar sector of the GSW model for g = 0, g' = gz = 2q, i.e. for sin^ 0^ = 1, 



and IF" = 0. 



Let us take a close look at the symmetries. The action is invariant under G = SU (2) global x 
U{l)iocai, with transformations 




$ ^ e*97(-)<i, ^1 11 ] ^ Y^^Y^ + a^7(x) , (52) 

under U{l)iocai, and 



^cos(^) +m3sin(^) 1(711 — m2)sin(^ 
1(711 +7?i2)sin(^) cos(^) — i7i3sin(^ 




under SU (2) global, where a = y^3?f + a| + af G [0, 47r) is a positive constant and ria = aa/a 
is a constant unit vector. Note that a shift of the function 7(2;) by 2Tr/q leaves the trans- 
formations unaffected. The model actually has symmetry G = [SU{2)giobai x U{^) local] /Z2; 
the Z2 identification comes because the transformation with (a, 7) is identified with that 
with (a + 27r, 7 + ir/q). Once $ acquires a vacuum expectation value, the symmetry breaks 
down to H = U{1) exactly as in the GSW model, except for the fact that the unbroken 
C/(l) subgroup is now global (for instance, if the VEV of the Higgs is ($"^) = 77(0, l)/^/2, 
the unbroken global U{1) is the subgroup with iii = 712 = 0, 713 = 1, (77 = a/2). Thus, the 
symmetry breaking is [SU{2)giobai x U{l)iocai]/Z2 U{l)giobai. 

Note also that, for any fixed $0 a global phase change can be achieved with either a global 
C^(l)iocai transformation or a >S'[/(2)giobai transformation. The significance of this fact will 
become apparent in a moment 

Like in the GSW model, the vacuum manifold is the three sphere 

2 

V = {$ eC^ I $t$ = !L} 53 ^ (54) 

which is simply connected, so there are no topological string solutions. On the other hand, if 
we only look at the gauged part of the symmetry, the breaking looks like C/(l) — > 1, identical 
to that of the Abelian Higgs model, and this suggests that we should have local strings. 

After symmetry breaking, the particle content is two Goldstone bosons, one scalar of mass 
nis = V2\r] and a massive vector boson of mass my — qrj.ln this section it will be convenient 
to use rescaled units throughout; after the rescaling (Isl), and dropping hats, we find 



q^S = / d^x 



\i^,^^Y,)^' - W.^Y^'^ - ^(<f>t$-l)' 



(55) 



and, as in the Abelian Higgs case, (3 = rri^/mf, = 2X/q^ is the only free parameter in the 
model. The equations of motion 

Df,D^'<i> + I3{\^\^ - 1)$ = 

(56) 

d^Yf,, ^ ~i<i>^ D.^ . 

are exactly the same as in the Abelian Higgs model but replacing the scalar field by the 
SU(2) doublet, and complex conjugation by hermitian conjugation of $. Therefore, any so- 
lution ^{x), Yfj_{x) of (^) (in rescaled units) extends trivially to a solution ^si{x), {Y^)si{x) 
of the semilocal model if we take 

<i>,i{x) = $(x)$o {Yf,)siix) = %{x) (57) 
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with $0 a constant SU{2) doublet of unit norm, ^J^o = 1. In particular, the Nielsen-Olesen 
string can be embedded in the semilocal model in this way. The configuration 

$ = /wo(p)e"'^$o, Y^nvNoip)dip (58) 

remains a solution of the semilocal model with winding number n provided /atq and vno 
are the solutions to the Nielsen-Olesen equations (^). In this context, the constant doublet 
$0 is sometimes called the 'colour' of the string (do not confuse with SU{3) colour!). One 
important difference with the Abelian Higgs model is that a scalar perturbation can remove 
the zero of $ at the center of the string, thereby reducing the potential energy stored in the 
core. 



Consider the energy per unit length, in these units, of a static, cylindrically symmetric 
configuration along the z-axis: 



E 



= / d'x 



i(5™r„ - dnY,nf + \{d„r - iF^O^P + ^($t$ » if 



(59) 



Note, first of all, that any finite energy configuration must satisfy 

{dm - iYm)(t>l ^ , (dm - iYm)(l>2 ^ , 01 01 + 02 </'2 ^ 1 aS p -> OO 

(As before, m,n = 1,2 and (p, <^) are polar coordinates on the plane orthogonal to the 
string). This leaves the phases of 0i and 02 undetermined at infinity and there can be 
solutions where both phases change by integer multiples of 2tt as we go around the string; 
however, there is only one U{1) gauge field available to compensate the gradients of 0i and 
02, and this introduces a correlation between the winding in both components: the condition 
of finite energy requires that the phases of 0i and 02 differ by, at most, a constant, as p ^ oo. 
Therefore, a finite energy string must tend asymptotically to a maximal circle on 



Y ndip 



or Y 



P 



(60) 



where < a < 1 and C are real constants, and determine the 'colour' of the string. A few 
comments are needed at this point. 



• Note that the choice of $o is arbitrary for an isolated string (any value of $o can be 
rotated into any other without any cost in energy) but the relative 'colour' between two 
or more strings is fixed. That is, the relative value of $o is significant whereas the absolute 
value is not. 



• The number n is the winding number of the string and, although it is not a topological 
invariant in the usual sense (the vacuum manifold, S'^, is simply connected), it is topolog- 
ically conserved. The reason is that, even though any maximal circle can be continuously 
contracted to a point on S^, all the intermediate configurations have infinite energy. The 
space that labels finite energy configurations is not the vacuum manifold but, rather, the 
gauge orbit from any reference point $0 G V, and this space [Giocai/ Hiocai), is not simply 
connected: it i{G local / H local) = ''^i{U{\)/\) — Thus, configurations with different wind- 
ing numbers are separated by infinite energy barriers, but this information is not contained 
in 7ri(V) 0. 

^The fact that the gauge orbits sit inside V = G/H without giving rise to non- 
contractible loops can be traced back to the previous remark that every point in 
the gauge orbit of can also be reached from with a global transformation. 
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• On the other hand, because tti{V) — 1, the existence of a topologicaUy conserved winding 
number does not guarantee that winding configurations are non-dissipative either. In con- 
trast with the Abehan Higgs model, a field configuration with non-trivial winding number 
at p = oo can be extended inwards for all p without ever leaving the vacuum manifold. 
Thus, the fact that TTi{Giocai / Hi^cai) 7^ 1 only means that finite energy field configurations 
fall into inequivalent sectors, but it says nothing about the existence of stable solutions 
within these sectors. 

• Thus, we have a situation where 

ni{V)^TriiG/H) = ni{S'') = l but ^i(Giocai/i?iocai) = ^1(5^) - £ • (61) 

and the effect of the global symmetry is to eliminate the topological reason for the existence 
of the strings. Notice that this subtlety does not usually arise because these two spaces are 
the same in theories where all symmetries are gauged (like GSW, Abelian Higgs, etc.). 
We will now show that, in the semilocal model, the stability of the string depends on the 
dynamics and is controlled by the value of the parameter /3 = 2X/q^. Heuristically we 
expect large /3 to mimic the situation with only global symmetries (where the strings would 
be unstable) , whereas small /3 resembles the situation with only gauge symmetries (where 
we expect stable strings). 



3.2 Stability 



Let us first prove that there are classically stable strings in this model. We can show this 



analytically for /3 = 1 |121 . Recall the expression of the energy per unit length (|59|). The 
analysis in the previous section goes through when the complex field is replaced by the 
SU{2) doublet, and we can rewrite 



E 

2TT\n\ 



j (fx |Di$±iZ?2*P + i(B±($t$_i))2 + i(^„i)($t$_i)2 ^ (62) 



choosing the upper or lower signs depending on the sign of n. Since n is fixed for finite 
energy configurations this shows that, at least for /3 = 1, a configuration satisfying the 
Bogomolnyi equations 

{Di±iD2)'^ = S±($^$-1)=0, (63) 

is a local minimum of the energy and, therefore, automatically stable to infinitesimal per- 
turbations. But these are the same equations as in the Abelian Higgs model, therefore the 
semilocal string (^ automatically saturates the Bogomolnyi bound (for any 'colour' $o)- 
Thus, it is classically stable for (3 = 1. 

This argument does not preclude zero modes or other configurations degenerate in energy. 
Hi ndmar sh ||55|| showed that, for /3 = 1 there are indeed such zero modes, described below 

dH). 



m 



We have just proved that, for /3 = 1, semilocal strings are stable. This is surprising because 
the vacuum manifold is simply connected and a field configuration that winds at infinity may 
unwind without any cost in potential enevgy^. The catch is that, because 'iTi{Giocai / Hiocai) — 

® In the Nielsen-Olesen case a configuration with a non-trivial winding number 
must go through zero somewhere for the field to be continuous. But here, a configura- 
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7ri([/(l)) = Z is non trivial, leaving the U{1) gauge orbit is still expensive in terms of 
gradient energy. 

As we come in from infinity, the field has to choose between unwinding or forming a semilocal 
string, that is, between acquiring mostly gradient or mostly potential energy. The choice 
depends on the relative strength of these terms in the action, which is governed by the value 
of and we expect the field to unwind for large (3, when the reduction in potential energy 
for going off the vacuum manifold is high compared to the cost in gradient energy for going 
off the U{1) orbits. And vice versa. Indeed, we will now show that, for (3 > 1, the n ~ 1 
vortex is unstable to perturbations in the direction orthogonal to $0 (5^] while, for /3 < 1, 
it is stable. For /3 = 1, some of the perturbed configurations become degenerate in energy 
with the semilocal vortex and this gives a (complex) one-parameter family of solutions with 
the same energy and varying core radius [q5| . 



3.2.1 The stability of strings with /3 > 1 

Hindmarsh has shown that for /3 > 1 the semilocal string configuration with unit 
winding is unstable to perturbations orthogonal to $0i which make the magnetic flux spread 



to infinity . As pointed out by Preskill [102|, this is remarkable because the total amount of 
flux measured at infinity remains quantized, but the flux is not confined to a core of finite 
size (which we would have expected to be of the order of the inverse vector mass). 

The semilocal string solution with n = 1 is, in rescaled units, 

= /ivo(p)e"^*o , Ysi = VNoip)dip . (64) 



However, as pointed out in [ p5| , this is not the most general static one- vortex ansatz com- 
patible with cylindrical symmetry. Consider the ansatz 

$ = /(p)e*'^$o+.9(p)e'"'^$± , Y = v{p)d^ , (65) 

with |$o| = I'J'il = 1 and (f>o$± = 0. The orthogonality of $0 and $^ ensures that the effect 
of a rotation can be removed from $ by a suitable SU{2) x U{1) transformation, therefore 
the configuration is cylindrically symmetric. For the configuration to have finite energy we 
require the boundary conditions /(O) = ^'(0) = w(0) = and / ^ 1, g ^ 0, v ^ 1 as 

p — !■ 00 

We know that if g = the energy is minimised by the semilocal string configuration / = 
/no J V — vjsfo, because the problem is then identical to the Abelian Higgs case. The question 
is whether a non-zero g can lower the energy even further, in which case the semilocal string 
would be unstable. The standard way to find out is to consider a small perturbation of ( |64[ ) 
of the form g — 0(p)e*"* and look for solutions of the equations of motion where g grows 
exponentially, that is, where < 0. The problem reduces to finding the negative eigenvalue 
solutions to the Schrodinger-type equation 

,p{p) = u^^Pip) (66) 

tion like $-^(/9 = oo) = r/(0, e*'^)/A/2 can gradually change to $'^(/9 = 0) = r/(l,0)\/2 
as we move towards the centre of the "string" without ever leaving the vacuum 
manifold. This is usually called 'unwinding' or 'escaping in the third dimension' by 
analogy with condensed matter systems like nematic liquid crystals. 
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Figure 2: A two-dimensional simulation of the evolution of a perturbed isolated semilocal 
string with /3 > 1, from [Q. The plot shows the (rescaled) energy density per unit length 
in the plane perpendicular to the string. /? = 1.1 The initial conditions include a large 
destabilizing perturbation in the core, = 0) = (1, fNo{p)e^'^), which is seen to destroy 
the string. 



Achucarro, Kuijken, Perivolaropoulos, Vachaspati 
Nucl. Phys. B388 (1992) 435, Fig. 4 



First of all, it turns out that it is sufhcient to examine the m ~ case only. Note that, 
since < v{p) < 1, for m > 1 the second term is everywhere larger than for m = 1, so if 
we one can show that all eigenvalues are positive for m — 1 then so are the eigenvalues for 
m > 1. But for m = 1 the problem is identical to the analogous one for instabilities in / 
in the Abelian Higgs model, and we know there are no instabilities in that case. Therefore 
it is sufficient to check the stability of the solution to perturbations with m — (negative 
values of m also give higher eigenvalues than m = 0.) 

If m = 0, the above ansatz yields 



E 



= 27r 



P 



iff + {g'f + 



l\2 



1 

2p 



(1 



9' 1) 



dp 
\67) 

for the (rescaled) energy functional (|59|). Notice that a non-zero g at p = (where / ^ 1) 
reduces the potential energy but increases the gradient energy for small values of p. If P is 
large, this can be energetically favourable (conversely, for very small /3, the cost in gradient 
energy due to a non-zero g could outweigh any reduction in potential energy). Indeed, 
Hindmarsh showed that there are no minimum-energy vortices of finite core radius when 
/3 > 1 by constructing a one-parameter family of configurations whose energy tends to the 
Bogomolnyi bound as the parameter po is increased: 



m = ^ 

Po 



1+^ 

Po. 



-1/2 



9{P) 



1 + ^ 
Po. 



-1/2 2 

Po 



1 + ^ 

Po 



(68) 



The energy per unit length of these configurations is E — Tirf'il -t- 1/3/9q) which, as po ^ oo, 
tends to the Bogomolnyi bound. This shows that any stable solution must saturate the 
Bogomolnyi bound, but this is impossible because, when /3 > 1, saturation would require 
B ~ Q everywhere, which is incompatible with the total magnetic flux being 27r/g (see 
the comment after eq. (|5^). While this does not preclude the possibility of a metastable 
solution, numerical studies have found no evidence for it |^. All indications are that, for 
/? > 1, the semilocal string is unstable towards developing a condensate in its core which 
then spreads to infinity. 



Thus, the semilocal model with /3 > 1 is a system where magnetic flux is quantized, the 
vector boson is massive and yet there is no confinement of magnetic fiux ^. 

^Preskill has emphasized that the "mixing" of global and local generators is a 



23 



Figure 3: The evolution of a string with f3 < 1. The initial configuration is the same as 
in Fig. H but now, after a few oscillations, the configuration relaxes into a semilocal string, 
= (0,/ivo(p)e^'^). /3 = 0.9 



Achucarro, Kuijken, Perivolaropoulos, Vachaspati 
Nucl. Phys. B388 (1992) 435, Fig. 1 



3.2.2 The stability of strings with /3 < 1 

Semilocal strings with /3 < 1 are stable to small perturbations. Numerical analysis of the 
eigenvalue equations p5| , p6| shows no negative eigenvalues, and numerical simulations of 
the solutions themselves indicate that they are stable to z-independent perturbations @, ^ , 
including those with angular dependence. Note that the stability to z-dependent perturba- 
tions is automatic, as they necessarily have higher energy. These results are confirmed by 
studies of electroweak string stability |]53|, 0| taken in the limit 6^ — > tt/2. 



3.2.3 (3 = 1 zero modes and skyrmions 

Substituting the ansatz ( |65| ) into the (rescaled) Bogomolnyi equations for n — 1 gives : 



P 

v'{p)+p{f{p)+9\p)-l) = Q 



When /3 = 1 we showed earlier that the semilocal string / — /nOi 9 — — vnq saturates 
the Bogomolnyi bound, so it is necessarily stable (since it is a minimum of the energy). 
There may exist, however, other solutions satisfying the same boundary conditions and 
with the same energy. Hindmarsh showed that this is indeed the case by noticing that the 
eigenvalue equation has a zero-eigenvalue solution pq| 



1^ = 1^0 exp 



P 



ipo ~ const , (70) 



which signals a degeneracy in the solutions to the Bogomolnyi equations. (Note that the 
'colour' at infinity, $0, is fixed, so this is not a zero mode associated with the global SU{2) 
transformations; its dynamics have been studied in fP^ .) 



necessary condition for this behaviour, that is, there must be a generator of H which 
is a non-trivial linear combination of generators of Ggiobal and Giocal [ |102|| . 
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It can be shown that the zero mode exists for any value oi g, not just 5 = 0; the Bogomolnyi 
equations (|6|) are not independent since, 

9{p) = (71) 
P 

is a solution of the second ec^uation for any (complex) constant qq. Solving the other two 
equations leads to the most general solution with winding number one and centred at p = 0. 
It is labelled by the complex parameter qq, which fixes the size and orientation of the vortex: 



^1 



^ ^ ^exp(iz.(p;|go|)^ (72) 



where m = In |<I>p is the solution to 

V2u + 2(l-e") = V2ln(p2 + |go|2)^ y as p ^ oo . (73) 



If qo 7^ 0, the asymptotic behaviour of these solutions is very different from that of the 
Nielsen-Olesen vortex; the Higgs field is non-zero at p = and approaches its asymptotic 
values like 0{p~^). Moreover, the magnetic field tends to zero as B 2\qo\'^p~'^, so the width 
of the flux tube is not as well-defined as in the qo = case when B falls off exponentially. 
These go 7^ solutions have been dubbed 'skyrmions'. In the limit \qo\ — > 0, one recovers 
the semilocal string solution (^), with u — ln(/^Q), the Higgs vanishing at p = and 
approaching the vacuum exponentially fast. On the other hand, when \qo\ ^ 1, u ~ the 
scalar field is in vacuum everywhere and the solution approximates a (CP^ lump js^, . 
Thus, in some sense, the 'skyrmions' interpolate between vortices and (DP^ lumps. 



3.2.4 Skyrmion dynamics 



We have just seen that, for /? = 1 the semilocal vortex configuration is degenerate in energy 
with a whole family of configurations where the magnetic fiux is spread over an arbitrarily 
large area. It is interesting to consider the dynamics of these 'skyrmions' when /3 7^ 1 |l^ : 
large skyrmions tend to contract if /? < 1 and to expand if /? > 1. The timescale for the 
collapse of a large skyrmion increases quadratically with its size . Thus large skyrmions 
collapse very slowly. 

Benson and Bucher p7| derived the energy spectrum of delocalized 'skyrmion' configurations 
in 2+1 dimensions as a function of their size. More precisely, they defined an 'antisize' 
X — -Emagnctic/^-totai as the ratio of the magnetic energy / (PxB^/2 to the total energy 
(pOf). Note that when the flux lines are concentrated, magnetic energy is high compared 
to the other contributions, and vice versa. Thus, x — > corresponds to the limit in which 
the magnetic flux lines are spread over an infinitely large area, which explains the name 
'antisize'. 



For large skyrmions - those with x < /3/(l + /3) - they concluded that the minimum energy 
configuration among all delocalized configurations with antisize x satisfies 

(if X > /3/(l + /3) the analysis does not apply). Therefore, energy decreases monotoni- 
cally with decreasing x for /? > 1 and increases monotonically for /3 < 1, confirming that 
delocalized configurations tend to grow in size if /3 > 1 and shrink if /3 < 1. 
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This behaviour is observed in numerical simulations Benson and Bucher have 
pointed out that in a cosmological setting the expansion of the Universe could drag the 
large skyrmions along with it and stop their collapse. The simulations in flat space are at 
least consistent with this, in that they show that delocalised configurations tend to live 
longer when artificial viscosity is increased, but a full numerical simulation of the evolution 
of semilocal string networks has not yet been performed and is possibly the only way to 
answer these questions reliably. 

Finally, we stress that the magnetic flux of a skyrmion does not change when it expands 
or contracts (the winding number is conserved) but this does not say anything about how 
localized the flux is. In contrast with the Abelian Higgs case, the size of a skyrmion can be 
made arbitrarily large with a finite amount of energy. 



3.3 Semilocal string interactions 



3.3.1 Multivortex solutions, (3 — 1, same colour 

Multi- vortex solutions in 2+1 dimensions corresponding to parallel semilocal strings with 
the same colour have been constructed by Gibbons, Ruiz-Ruiz, Ortiz and Samols [ ^Sf for 
the critical case (3=1. Their analysis closely follows that of |64| in the case of the Abelian 
Higgs model, and starts by showing that, as in that case, the full set of solutions to the 
(second order) equations of motion can be obtained by analysing the solutions to the (first 
order) Bogomolnyi equations. 

In the Abelian Higgs model, solutions with winding number n are labelled by n unordered 
points on the plane (those where the scalar field vanishes) which, for large separations, are 
identified with the positions of the vortices. In the semilocal model, the solutions have other 
degrees of freedom, besides position, describing their size and orientation. 

Assuming without loss of generality that the winding number n is positive, and working in 
temporal gauge = 0, any solution with winding number n is specified (up to symmetry 
transformations) by two holomorphic polynomials 

n 
r=l 

= z" + Pn_iz"^^ + . . . + piZ + po 

and 

Q„(z) = qn-iz''-^ + . . . + qiz + qo (75) 

where z = x + iy is a complex coordinate on the xy plane. The solution for the Higgs fields 
is, up to gauge transformations, 

f(j)i\_ (ih<^^) 

where the function u{z,'z) = ln(|0ip + |(/)2p) must satisfy 

V^u + 2(1 - e") - ln(|P„|2 + |Q„|2) , (77) 

and tend to as \z\ — > oo. Although its form is not known explicitly, ref. proved the 
existence of a unique solution to this equation for every choice of P„ and Q„ (if P„ and Qn 
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Figure 4: A numerical simulation of the interaction between two parallel semilocal strings 
with different 'colour', from Ref. M. The initial configuration has one string with = 
(0, /(pi)e*'^^) and the other with $2 = {if{p2)e^'^^,0), where {pi,ipi) are polar coordinates 
centred at the cores of each string. The energy density of the string pair is plotted in the 
plane perpendicular to the strings. The colour difference is radiated away in the form of 
Goldstone bosons, and the strings cores remain at their initial positions. f3 = 0.5. 



have a common root then exp[u/2] has a zero there, so the expression for the Higgs field is 
everywhere well-defined) . The gauge field can then be read off from the Bogomolny equations 
(|63|). This generalises ( [72| ) to arbitrary n. The coefficients of Pn{z), Qn{z) parametrise the 
moduli space, C^". 

The Nielsen-Olesen vortex has Qr, 
finds 



0. If P„ 7^ 0, then in regions where \Qn\ << \Pn\ one 



f - 



Qn 



Pn 
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Pn 



1 - 



Qn 



Pn 



(78) 



indicating that the scalar fields fall off as a power law, as opposed to the usual exponential 
fall off found in NO vortices. The same is true of the magnetic field. 



The low energy scattering of semilocal vortices and skyrmions with (3—1 was studied in 
[ [tsI in the geodesic approximation of |Q . The behaviour of these solitons was found to be 
analogous to that of (CP^ lumps but without the singularities, which are smoothed out in 
the core. 



3.3.2 Interaction of parallel strings, (3 <1, different colours 

Ref. carried out a numerical study in two dimensions of the interaction between stable 
(/3 < 1) strings with different "colour" with non-overlapping cores. It was found that the 
strings tend to radiate away their colour difference in the form of Goldstone bosons, and 
there is little or no interaction observed. The position of the strings remains the same during 
the whole evolution while the fields tend to minimize the initial relative SU (2) phase (see 
figure H). 

Thus, we expect interactions betwen infinitely long semilocal strings with different colours 
to be essentially the same as for Nielsen-Olesen strings. This expectation is confirmed by nu- 
merical simulations of two- and three-dimensional semilocal string networks p, W] , discussed 



in B.f; 
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3-4 Dynamics of string ends 



Note that, in contrast with Nielsen-Olesen strings, there is no topological reason that forces 
a semilocal string to continue indefinitely or form a closed loop. Semilocal strings can end 
in a "cloud" of energy, which behaves like a global monopole |^ . 

Indeed, consider the following asymptotic configuration for the Higgs field: 



V2 V sini^ 

which is ill-defined &t 9 ~ it and at r = 0. We can make the configuration regular by 
introducing profile functions such that the Higgs field vanishes at those points: 

?/ / /ii(r, 6*) cos i( 



where hi and /i2 vanish at r = and h2{r, tt) = 0. This configuration describes a string in 
the z < axis ending in a monopole at z = 0. 

At large distances, r >> 1, the Higgs field is everywhere in vacuum (except at w tt) 
and we find <i>V$ ~ x, just like for a Hedgehog in 0(3) models. On the other hand, 
the configuration for the gauge fields resembles that of a semi-infinite solenoid; the string 
supplies U(l) flux which spreads out from z = Q. 

This is the 6^ ^ tt/2 limit of a configuration first discussed by Nambu |Q in the context 
of the GSW model -see section |^ - but here the energy of the monopole is linearly divergent 
because there are not enough gauge fields to cancel the angular gradients of the scalar field. 

Angular gradients provide an important clue to understand the dynamics of string ends. If 
/3 < 1, numerical simulations show that string segments grow to join nearby segments or 
to form loops (see figures || and ^) This confirms analytical estimates in refs. [|8[ 
In other cases the string segment collapses under its own tension, with the monopole and 
antimonopole at the ends annihilating each other. 



3.5 Numerical simulations of semilocal string networks 



As the early Universe expanded and cooled to become what we know today it is very 
likely that it went through a number of phase transitions where topological (and possibly 
non-topological) defects are expected to have formed according to the Kibble mechanism 
Jn] , 134, |105[ |. Although the cosmological evidence for the existence of such defects remains 
unclear M, there is plenty of experimental evidence from condensed matter systems that 
networks of defects do form in symmetry breaking phase transitions [ p4[ , the most recent 
confirmation coming from the Lancaster-Grenoble-Helsinki experiments in vortex formation 
in superfluid Helium [ p5[ . An important question is whether semilocal (and electroweak) 
strings are stable enough to form in a phase transition. 



We defer discussion of the electroweak case to section \)A\ Here we want to review recent 
numerical simulations of the formation and evolution of a network of /? < 1 semilocal strings 
^, 1^ which show that such strings should indeed form in appreciable numbers in a phase 
transition. The results suggest that, even if no vortices are formed immediately after <& 
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Figure 5: Loop formation from semilocal string segments. The figure shows two snapshots, 
at i = 70 and t = 80, of a 64'^ numerical simulation of a network of semilocal strings with 
P = 0.05 from Ref_H, where the ends of an open segment of string join up to form a closed 
loop (see section |3.5| for a discussion of the simulations) . Subsequently the loops seem to 
behave like those of topological cosmic string, contracting and disappearing. 



has acquired a non-zero vacuum expectation value, the interaction between the gauge fields 
and the scalar fields is such that vortex formation does eventually occur simply because it 
is energetically favourable for the random distribution of magnetic fields present after the 
phase transition to become concentrated in regions where the Higgs field has a value close 
to that of the symmetric phase. 

Even though they do not account for the expansion of the Universe, these simulations 
represent a first step towards understanding semilocal string formation in cosmological 
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Figure 6: The growth of string segments to form longer strings. The figure shows two 
snapshots, at time i = 60 and i = 70 of a large 256'^ numerical simulation of a network 
of semilocal strings with [3 = 0.05 from Ref. Note several joinings of string segments, 
e.g. two separate joinings on the long central string, and the disappearance of some loops. 
The different apparent thickness of strings is entirely an effect of perspective. The simulation 
was performed on the Cray T3E at the National Energy Research Scientific Computing 
Center (NERSC). See section 3.5 for a discussion. 



phase transitions and they have already provided very interesting insights into the dynamical 
evolution of such a network. 
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3.5.1 Description of the simulations 



^From a technical point of view, the numerical simulation of a network of semilocal strings 
has additional complications over that of U{1) topological strings. Because there are not 
enough gauge degrees of freedom to cancel all of the scalar field gradients, the existence of 
string cores depends crucially on the way the fields (scalar and gauge) interact. Another 
problem, generic to all non-topological strings, is that the winding number is not well 
defined for configurations where the scalar is away from a maximal circle in the vacuum 
manifold, and this makes the identification of strings much more difficult than in the case 
of topological strings. 

The strategy proposed in to circumvent these problems is to follow the evolution of 
the gauge field strength in numerical simulations, since the field strength provides a gauge 
invariant indicator for the presence of vortices. The initial conditions are obtained by an 



extension of the Vachaspati-Vilenkin algorithm [124| appropriate to non-topological defects, 
plus a short period of dynamical evolution including a dissipation term (numerical viscosity) 
to aid the relaxation of configurations in the 'basin of attraction' of the semilocal string. 

As with any new algorithm, it is essential to check that it reproduces previously known 
results accurately, and this has been done in Note that setting 02 in the semilo- 
cal model obtains the Abelian Higgs model, thus comparison with topological strings is 
straightforward, and it is used repeatedly as a test case, both to check the simulation tech- 
niques and to minimise systematic errors when quoting formation rates. In particular, the 
proposed technique is tested in a two-dimensional toy model (representing parallel strings) 
in three different ways: a) restriction to the Abelian Higgs model gives good agreement with 



analytic and numerical estimates for cosmic strings in |124|; b) the results are robust under 
varying initial conditions and numerical viscosities (see Figure ^ , and c) they are in good 
agreement with previous analytic and numerical estimates for semilocal string formation in 

The results are summarized in Fig. |^. We refer the reader to refs. ^, |^ for details; 
however, a few comments are needed to understand those figures. 



• The study takes place in flat spacetime. Temporal gauge and rescaled units (^2|) are 
chosen. Gauss' law, which here is a constraint derived from the gauge choice Iq = 0, is used 
to test the stability of the code. 

• Space is discretized into a lattice with periodic boundary conditions. The equations of 
motion (^ ) are solved numerically using a standard staggered leapfrog method; however, to 
reduce its relaxation time an ad hoc dissipation term was added to each equation (77$ and 
rfYi respectively). A range of strengths of dissipation was tested, and it did not significantly 
affect the number densities obtained. The simulations displayed in this section all have have 
77 = 0.5. 



The number density of defects is estimated by an extension of the Vachaspati-Vilenkin 



algorithm [124| by first generating a random initial configuration for the scalar fields drawn 
from the vacuum manifold, which is not discretised, and then finding the gauge field con- 
figuration that minimizes the energy associated with (covariant) gradients^. If space is a 
grid of dimension , the correlation length is chosen to be some number p of grid points 
(p = 16 in P, B; the size of the lattice is either iV = 64 or = 256.) To obtain a reasonably 



lOi 



'^In fact, it turns out that the energy-minimization condition is redundant, since 
the early stages of dynamical evolution carry out this role anyway. 
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smooth configuration for the scalar fields, one throws down random vacuum values on a 
(N/p)^ subgrid; the scalar field is then interpolated onto the full grid by bisection. Strings 
are always identified with the location of magnetic flux tubes. 

For cosmic strings, the two-dimensional toy model accurately reproduces the formation rates 
of |124[ . For semilocal strings, on the other hand, the initial configurations generated in this 
way have a complicated flux structure with extrema of different values (top panel of Fig. 

, and it is far from clear which of these, if any, might evolve to form semilocal vortices; 
in order to resolve this ambiguity, the initial configurations are evolved forward in time. 
As anticipated, in the unstable regime P > 1 the fiux quickly dissipates leaving no strings. 
By contrast, in the stable regime P < 1 stringlike features emerge when configurations in 
the "basin of attraction" of the semilocal string relax unambiguously into vortices (bottom 
panel of Fig. |). 

Since the initial conditions are somewhat artificial, the results were checked against various 
other choices of initial conditions, in particular different initial conditions for the gauge field 
and also thermal initial conditions for the scalar field (see Fig. ||). All the initial conditions in 
] had zero initial velocities for the fields. Initial conditions with non-zero field momenta 
have not yet been investigated. 



3.5.2 Results and discussion 

These simulations give very important information on the dynamics and evolution of a 
network of semilocal strings. In particular, they confirm our discussion in the previous 
subsection of the behaviour of the ends of string segments, and of strings with different 
colours. String segments are seen to grow in order to join nearby ones or form closed loops, 
and very short segments are also observed to collapse and disappear. The colour degrees of 
freedom do not seem to introduce any new forces between strings. Because the strings tend 
to grow or form closed loops, time evolution makes the network resemble more and more a 
network of topological strings (NO vortices) but with lower number densitie^^ 

Note that the correlation length in the simulations is constrained to be larger than the size 
of the vortex cores, to avoid overlaps. This results in a minimal value of the parameter /3 
of around 0.05 (if P is lowered further, the scalar string cores become too wide to fit into 
a correlation volume, in contradiction with the vacuum values assumed in a Vachaspati- 
Vilenkin algorithm). Figure ^ shows the results for seven different values of /3 by taking 
several initial configurations on a 64'^ grid smoothed over every 16 grid-points. As expected, 
for /3 < 1 the formation rate depends on /3, tending to zero as P tends to 1. The ratio of 
semilocal string density to cosmic string density in an Abelian Higgs model for the same 
value of P is less than but of order one. For the lowest value of P simulated (/3 = 0.05), the 
semilocal string density is about one third of that of cosmic strings. 

One final word of caution about the possible cosmological implications of these simulations. 
We mentioned above that numerical viscosity was introduced to aid the relaxation of con- 
figurations close to the semilocal string. In an expanding Universe the expansion rate would 
provide some viscosity, though rj would typically not be constant. This may have an impor- 
tant effect on the production of strings. Indeed, note the different numbers of upward and 
downward pointing flux tubes in Fig. ^ despite the zero net flux boundary condition. The 
missing flux resides in the smaller 'nodules', made long-lived by the numerical viscosity; 

^^However, one important point is that no intersection events were observed in the 
semilocal string simulations, so the rate of reconnection has not been determined. 
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Figure 7: The flux tube structure in a two-dimensional semilocal string simulation with 

f3 = 0.05, from Ref ||^. The upper panel {t = 0) shows the initial condition after the process 

described in the text. The lower panel shows the configuration resolved into five flux tubes 

by a short period of dynamical evolution {t = 100). These flux tubes are semilocal vortices. 
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Figure 8: A test of the sensitivity of the results to the choice of initial conditions in 



a two-dimensional simulation with the algorithm proposed in section 3.5, The plot shows 



the number of semilocal strings formed per initial two-dimensional correlation volume. Each 
point is an average over ten simulations. Squares indicate that the vacuum initial conditions 
described in the text were used, while open circles indicate that non-vacuum (thermal) 
initial conditions were used. Both sets of initial conditions are seen to give comparable 
results. Statistical results are derived from a large suite of simulations (700 in all) carried 
out on a 64'^ grid (from Ref. 0) 



these are none other than the 'skyrmions' described in section As was explained there, 
the natural tendency of skyrmions when /3 < 1 is to collapse into strings, but the timescale 
for collapse increases quadratically with their size and Benson and Bucher have argued 
that the effect of the expansion could stop the collapse of large skyrmions almost completely. 
On the other hand one expects skyrmions to be formed with all possible sizes, so the ef- 
fect of the expansion on the number density of strings remains an open question. Another 
important issue that has not yet been addressed is whether these semilocal networks show 
scaling behaviour, and whether reconnections are as rare as the above simulations suggest. 
Both would have important implications for cosmology. However, the answer to these and 
other questions may have to wait until full numerical simulations are available. 



3.6 Generalisations and final comments 



i) Charged semilocal vortices 

The semilocal string solution described earlier in this section is strongly static and z- 
independent, by which we mean that Dt($) = i?z(<&) — 0. It is possible to relax these 
conditions while still keeping the Lagrangian and the energy independent of z The idea 
is that, as we move along the z-direction, the fields move along the orbit of the global 
symmetries; in other words, Goldstone bosons are excited. 

Abraham has shown that it is possible to construct semilocal vortices with finite energy 
per unit length carrying a global charge [|l| in the Bogomolnyi limit /3 = 1 They satisfy 

"'^^By contrast, charged solutions with Dq{^) 7^ in the Abelian Higgs model have 
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Figure 9: The ratio of lengths of semilocal and cosmic strings as a function of the stability 
parameter /3, from B. 



a Bogomolnyi-type bound and are therefore stable. Perivolaropoulos |100| has constructed 
spinning vortices (however these have infinite energy per unit length) . 

ii) Semilocal models with S'[/(A^)giobai x C^(l)iocai symmetry 

The generalization of semilocal strings to so-called Extended Abelian Higgs model s wit h an 



N-component multiplet of scalars whose overall phase is gauged is straightforward 1 121 , ^ , 
and has been analysed in detail in j5^, |4^. The strings are stable (unstable) for /3 < 1 
{f3 > 1) and for /3 = 1 they are degenerate in energy with skyrmionic configurations labelled 
by an iV — 1 complex vector. For winding n, and widely separated vortices, the Nn complex 
parameters that characterize the configurations can be thought of as the n positions in 
~(C and the {N - l)n 'orientations'. 

iii) Semilocal monopoles and generalized semilocality 

We have seen that semilocal strings have very special properties arising from the fact that 
ni{G/H) = but 7ri(Giocai/ffiocai) 7^ 0. An immediate question is whether it is possible to 
construct other non-topological defects such that 

7rfc(G/i/)=0 but 7rfc(Glocal/iflocal) ^ . (81) 

This possibility would be particularly interesting in the case of monopoles, k — 2, since they 
might retain some of the features of global monopoles, in particular a higher annihilation rate 
in the early Universe. Surprisingly, the answer seems to be negative. Within a very natural 



set of assumptions, it was shown in [121| that the condition ( |8l| ) can only be satisfied if the 
gauge group Giocai is Abelian, and therefore one cannot have semilocal monopoles (nor any 
other defects satisfying conditions ( ^l|) with fc > 1). 

However, Preskill has remarked that it is possible to define a wider concept of semilocality 



[102 1 by considering the larger approximate symmetry Gapprox which is obtained in the limit 
where gauge couplings are set to zero. The symmetry Gapprox is partially broken to the exact 
symmetry G ~ Giocai x G global (modulo discrete transformations) when the gauge couplings 



infinite energy per unit length |67] 
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are turned on It is then possible to have generahzed semilocal monopoles associated with 
non-contractible spheres in Giocai/-ff local which are contractible in the approximate vacuum 
manifold Gapprox/^^approx even though they are still non-contractible in the exact vacuum 
manifold G/H. 

Another obvious possibility is to have topological monopoles with "colour" , by which we 
mean extra global degrees of freedom, if the symmetry G ^ Ggiobai x Giocai is such that 
the gauge orbits are non-contractible two-spheres, 7r2(Giocai/-ffiocai) 7^ 1- Given that there 



are no semilocal monopoles |121], these monopoles must have 7r2 (G/iJ) 7^ 1, so they are 



topologically stable, and they have additional global degrees of freedom. 

iv) Semilocal defects and Hopf fibrations 

In the semilocal model, the action of the gauge group fibres the vacuum manifold as 
a non-trivial bundle over 5*2 ~ (CP\ the Hopf bundle. The fact that this bundle IS non- 
trivial is at the root of conditions (|6l|), and is ultimately the reason why the topological 
criterion for the existence of strings fails. In view of this, Hindmarsh fs^ has proposed an 
alternative definition of a semilocal defect: it is a defect in a theory whose vacuum manifold 
is a non-trivial bundle with fibre Giocai/-f^iocai- 

Extended Abelian Higgs models are similarly related to the fibrations of the odd- 
dimensional spheres S"^^^^ with fibre and base space (CP^~^. A natural question to ask 
is if the remaining Hopf fibrations of spheres can also be realised in a field theoretic model. 
This question was answered affirmatively in fs?! for the S"^ 5'^ fibration in a quaternionic 
model. Other non-trivial bundles were also implemented in this paper, but to date the field 
theory realisation of the S^^ ^ Hopf bundle remains an open problem. 

v) Monopoles and textures in the semilocal model: 

Since the gauge field is Abelian, divB = 0, and isolated magnetic monopoles are necessarily 
singular in semilocal models. The only way to make the singularity disappear is by embed- 
ding the theory in a larger non- Abelian theory which provides a regular core, or by putting 
the singularity behind an event horizon . One important question that has not yet been 
addressed is if the scalar gradients in these spherical monopoles make them unstable to 
angular collapse into a flux tube. A related system where this happens is in 0(3) global 
monopoles where the spherically symmetric configuration is unstable. In the semilocal case, 
it is possible that the pressure from the magnetic field might prevent the instability towards 
angular collapse. 

Finally, note that, because 7r3(S'"^) = E, there is also the possibility of textures in the 



semilocal model (51). In contrast with purely scalar 0(4) models, their collapse seems to 
be stopped by the pressure from the magnetic field . Of course they can still unwind by 
tunnelling. 

vi) We should point out that systems related to the semilocal model have been studied in 
condensed matter. In [E5|, the system was an unconventional superconductor where the role 



of the global SU(2) group was played by the spin rotation group. In |12J ] the hypothetical 
case of an "electrically charged" A-phase of '^He, i.e. a superconductor with the properties 



of '^He-A, was considered (see section 11.1 for a brief discussion of the A and B phases of 
■^He). In this case the global group was SO (3), the group of orbital rotations. Both papers 
discussed continuous vortices in such superconductors, which correspond to the "skyrmions" 
discussed here. 
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4 Electroweak strings 



In this section we introduce electroweak strings. There are two kinds: one, more precisely 
known as the Z-string, carries Z- magnetic flux, and is the type that was discussed by Nambu 
and that becomes stable as it approaches the semilocal limit. It is associated with the 
subgroup generated by 

Tz = n^T"" - sin^ O^Q . 

There are other strings in the GSW model that carry SU{2) magnetic flux, called W- 
strings. There is a one-parameter family of W strings which are all gauge equivalent to one 
another, and they are all unstable. They are generated by a linear combination of the SU(2) 
generators T'^ and T^. These will be discussed in more detail in the next section. 



4-1 The Z string 



Modulo gauge transformations, the configuration describing a straight, infinitely long Z- 



string along the z-axis is 117 



A„ 



Z = VNoip)dip or, in vector notation, Z = 

9z 







vncAp) 
P 



-LP 



(82) 



where / and v are the Nielsen-Olesen profiles that solve the equations (|4C|). It is straightfor- 
ward to show that this is a solution of the bosonic equations of motion (alternatively, one 
can show that it is an extremum of the energy [117|). Equations ( p2| ) describe a string with 
unit winding. The solutions with higher winding number can be constructed in an analogous 
way, but note that the winding number is not a topological invariant. The unstable string 
can decay by unwinding until it reaches the vacuum sector. 



The solution (^) reduces to the semilocal string in the limit sin^ 9^ — 1, and therefore it is 
classically stable for /3 < 1 and unstable for (3 > 1 (see section where (3 is now the ratio 
between the Higgs mass, \/2\r\ and the Z-boson mass gzVl'^^ thus 



8A 



(83) 



The Z-string configuration is axially symmetric, as it is invariant under the action of the 
generalised angular momentum operator 

^L, + S,+h (84) 
where L^, Sz and Iz are the orbital, spin and isospin parts, respectively, defined in section 



9.1 



The Z-string carries a Z-magnetic fiux 



Fz = — (85) 
9z 
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thus particles whose Z charge is not an integer multiple of 5z/2 will have Aharonov-Bohm 
interactions with the string (see section 8.3). The Z-string can terminate on magnetic 
monopoles (such configurations are discussed in section]^). When a string terminates, the 
discrete Aharonov-Bohm interaction can be smoothly deformed to the trivial interaction. 
The smoothness is provided by the presence of the magnetic flux of the monopole. 



• Note that, in the background given by (p^), the covariant derivative becomes 

= i?^|z-string - + *y [-'^{T^ ~ Q Sin' 0.^)] (86) 

in particular, left and right fermion fields couple to with different strengths, since the 
effective Z-charge 

ci=-2{T^ -Qsin^e^) (87) 

has different values, qR = ± 1. (Note that q is proportional to the string generator T^, 
defined in equation (p^. The proportionality factor has been introduced for later conve- 
nience). This will be important in section |8.3| . Note also that, for the Higgs field, 

q = diag(-cos26'^, 1) (88) 



• Ambj0rn and Olesen |10| and, more recently, Bimonte and Lozano [ pO[ have derived 
Bogomolnyi-type bounds for periodic configurations in the GSW model. They consider static 
configurations such that physical observables are periodic in the a;?/-plane and cylindrically 
symmetric in each cell. If A is the area of the basic cell, they find that the energy (per unit 
length) satisfies 

E > (l/2e2)m2,(2g'Fy - A) if mn > mz 

(89) 

E > {l/2e^){m.wmH/mzf{2g'FY - m^^A) if < mz 

where Fy is the magnetic flux of the hypercharge field through the cell. Note that the top 
line of ( |89| ) reduces to the familiar E > ($'l'$)gJV for the Abelian Higgs and semilocal case 
in the 5 — > limit (with q = g' /2). In the non- Abelian case the bound involves an area term 
and therefore does not admit a topological interpretation. 

In the Bogomolnyi limit, mn — mz, the bound is saturated for configurations satisfying 
the first order Bogomolnyi equations 

Di + iD2<S> = 

W12 + I^V''* = 

A solution to these equations describing a lattice of Z-strings was constructed in . Other 
periodic configurations with symmetry restoration had been previously found in the presence 
of an external magnetic field in pO| | . 



5 The zoo of electroweak defects 



The electroweak Z-string is one member in the zoo of electroweak defects. Other members 
include the electroweak monopole, dyon and the W-string. The latter fall in the class of 
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"embedded defects" and this viewpoint provides a simple way to characterize them. The 
electroweak sphaleron is also related to the electroweak defects. 



5.1 Electroweak monopoles 



To understand the existence of magnetic monopoles in the GSW model, recall the following 
sequence of facts: 

• The Z-string does not have a topological origin and hence it is possible for it to 
terminate. 

• As the hypercharge component of the Z-field in the string is divergenceless it cannot 
terminate. Therefore it must continue from within the string to beyond the terminus. 

• However, beyond the terminus, the Higgs is in its vacuum and the hypercharge mag- 
netic field is massive. Then, if the massive hypercharge flux was to continue beyond 
the string, it would cost an infinite amount of energy and this is not possible. 

• The only means by which the hypercharge field can continue beyond the terminus is 
in combination with the SU(2) fields such that it forms the massless electromagnetic 
magnetic field. 

So the terminus of the Z-string is the location of a source of electromagnetic magnetic field, 
that is, a magnetic monopole fo^ . We now make this argument more quantitative. 

Assume that we have a semi- infinite Z-string along the —z axis with terminus at the origin 
(see Fig. |l^). Let us denote the A- and Z- magnetic fluxes through a spatial surface by Fa 
and Fz ■ These are given in terms of the W- and Y- fluxes by taking surface integrals of the 
field strengths (see eqs. (|T9|), (|20|)). Therefore 

Fz = cosO^Fn - sine^FY , Fa = sin 6l^i^„ -f cos 6^,Fy , (91) 

where we have denoted the SU(2) flux (parallel to in group space) by F„ and the 
hypercharge flux by Fy. 

Now consider a large sphere E centered on the string terminus. The field configuration is 
such that there is only A-flux through S except near the South pole (S) of E, where there 
is only a Z magnetic flux. Hence, 

Fzb-s = , Fa\s = . (92) 

Together with (|9|) this gives, 

F„|s-s = tan0„Fy|E-s , F^\s ^- cot 9^Fy\s ■ (93) 

The hypercharge flux must be conserved as it is divergenceless. So 

Fy\s~s ^ -Fy\s = Fy , (94) 
and, inserting this and (^3|) in (|9|) yields 

J^a|s-s = -^, Fz\s^^^. (95) 

cos dyj Sm dyj 

Now the flux in the string along the —z axis is quantized in units of 'i'K/gz (recall 
Qz = e/ cos 6w sin6'u, gives the coupling of the Z boson to the Higgs fleld). Therefore, for the 
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Figure 10: The outgoing hypercharge flux of the monopole passing through the surface 
T, — S should equal the incoming hypercharge flux through the Z-string. 




unit winding string, 
Then (PI) yields, 



Fz\s^ — . (96) 
9z 



4:7T Att 47r o 

Fy = — sin6'u, , Fa\s-s = — tan6'^ = — sin 6^ (97) 

9z 9z e 

Hence the terminus of the string has net A-flux emanating from it and hence it is a magnetic 
monopole. 

The electromagnetic flux of the electroweak monopole appears to violate the Dirac quanti- 
zation condition. However this is not true since one must also take the Z-string into account 
when deriving the quantization condition relevant to the electroweak monopole This be- 
comes clearer when we work out the magnetic flux for the SU{2) fields. Using (|9^ ) with 
(p7|), the net non-Abelian flux is: 

47r 

Fn ^ Fn\s + Fnl^^s ^ — (98) 
9 

just as we would expect for a 't Hooft-Polyakov monopole That is, the Dirac quantiza- 
tion condition works perfectly well for the SU(2) field and the monopole charge is quantized 
in units of -Att / g. Another way of looking at ( |98| ) is to say that the electroweak monopole is a 
genuine SU(2) monopole in which there is a net emanating C/(l)„ C SU{2) flux. The struc- 
ture of the theory, however, only permits a linear combination of this flux and hypercharge 
flux to be long range and so there is a string attached to the monopole. But this string is 
made of Z field which is orthogonal to the electromagnetic field and so the string does not 
surreptitiously return the monopole electromagnetic flux. Also, the magnetic charge on the 
monopole is conserved and electroweak monopoles can only disappear by annihilating with 
antimonopoles . 

It is useful to have an explicit expression describing the asymptotic field of the electroweak 
monopole and string. Nambu's monopole-string configuration, denoted by ($, W^,y^), is 



^^^{ cos(0/2) 
V2 \sm{e/2)e^^ 



(99) 



where, 9 and (p are spherical coordinates centred on the monopole, and the gauge field 
configuration is, 

gW;^ = -e°^'=n^9^n'= -f- icos^ 6'u,n"($t _ a^$t |,) (iqo) 
g% = -i sin^ 0„($t _ a^$t $) (101) 



where, n" is given in eq. (|16|). 

Note that there is no electroweak configuration that represents a magnetic monopole sur- 
rounded by vacuum. 



5.2 Electroweak dyons 



Given that the electroweak monopole exists, it is natural to ask if dyonic configurations exist 
as w ell. We now write down dyonic configurations that solve the asymptotic field equations 
[120|. The existence of such configurations is implicit in Nambu's original paper in the guise 
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of what he called "external" potentials |Q . Essentially, the dyon solution is an electroweak 
monopole together with a particular external potential. 

The ansatz that describes an electroweak dyon connected by a semi- infinite Z string is: 

$ = # (102) 

- dt ^ (103) 

cos 

Y = Y - dt^— (104) 
sin 9m 

where, C = ({t,x), overdots denote partial time derivatives and barred fields have been 
defined in the previous subsection. 



We now need to insert this ansatz into the field equations and to find the equation satisfied 
by (. Some algebra leads to 

d'd^C = , dtdX = (105) 
which can be solved by separating variables, 



This leads to 



(106) 



(107) 



The particular solution that we will be interested in is the solution that gives a dyon. Hence, 
we take: 

c c ^ ff \ <? sin 6*^ cos 6*^ 1 

f = Coi , f{r) = —z , (108) 

47r^o r 



where, and q are constants. Now, using (108), together with ( |103| ), (104) and ( |106D , we 
get the dyon electric field: 

Ea = ^^. (109) 
An 



For a long segment of string, the monopole and the antimonopole at the ends are well 
separated and we can repeat the above analysis for both of them independently. Therefore, 
the electric charge on the antimonopole at one end of a Z-string segment is uncorrelated 
with the charge on the monopole at the other end of the string. This means that we can 
have dyons of arbitrary electric charge at either end of the string. The situation will change 
with the inclusion of fermions since these can carry currents along the string and transport 
charge from monopole to antimonopole. 

This completes our construction of the dyon-string system in the GSW model. As of now, 
the charge q on the dyon is arbitrary. Quantum mechanics implies that the electric charge 
must be quantized. If wc include a 9 term in the electroweak action (but no fermions): 

So-^l d'xW^^^W^''^ (110) 

where 

^M^" = IgM^^'-W-;^ , (111) 
then the charge quantization condition becomes 




(112) 
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This agrees with the standard result for dyons ]133| |. 



In the GSW model with fcrmions, it is known that the 9 term can be eliminated by a 
rotation of the fermionic fields. This argument can be turned around to argue that the 
CP violation in the mass matrix of the fermions will lead to an effective 9 term and so 
the electroweak monopoles will indeed have a fractional charge with 9 being related to 
the CP violation in the mass matrix. The precise value of the fractional electric charge on 
electroweak monopoles has not yet been calculated and remains an open problem. 

It should be mentioned that, even though the electric charge on an electroweak dyon can be 



fractional as in (112), the total electric charge on the dyon-string system is always integral 
because the CP violating fractional charge on the monopole is equal and opposite to that 
on the antimonopole. 



5.3 Embedded defects and W-strings: 



A very simple way of understanding the existence of electroweak string solutions is in terms 
of embedded defects. While this method does not shed any light on the stability of the 
electroweak string, it does provide a scheme for finding other solutions. 

The idea is that the electroweak symmetry group contains several U{1) subgroups which 
break completely when the electroweak symmetry breaks. Corresponding to each such break- 
ing, one m ight have a string solution. A more complete analysis tells us when such a solution 
can exist [p| || 0. 



Consider the general symmetry breaking 

G^H (113) 

Suppose Gemb IS a subgroup of G which, in this process, breaks down to Gemb H H. Then 
we ask the question; when are topological defects in the symmetry breaking 

Gemb Gemb H H (114) 

also solutions in the full theory? An answer to this question requires separating the gauge 
fields into those that transform within the Gemb subgroup and those that do not. Similarly, 
the Higgs field components are separated into those that lie in the embedded vector space 
of scalar fields and those that do not. Then, it is possible to write down general conditions 
under which solutions can be embedded @ ||. Here we shall not describe these conditions 
but remark that the Z-string is due to the embedded symmetry breaking 

Uil)z^l (115) 

where the U{l)z is generated by Tz, defined in eq. (p^). Now, there are other J7(l)'s that 
can be embedded in the GSW model which lie entirely in the SU (2) factor. For example, we 
can choose U{l)i which is generated by (one of the off-diagonal generators of SU{2)). 
Since we have 

C/(l)i ^ 1 (116) 

when the electroweak symmetry breaks, there is the possibility of another string solution in 
the GSW model. Indeed, it is easily checked that this string can be embedded in the GSW 
model and the solution is called a W-string. By considering a one parameter family of U{1) 
subgroups generated by 

= cos(C) + sin(C) (117) 
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we can generate a one parameter family of W-strings 



*=>o(.)(,,»X^) (118, 



2 2 

— — cos C V]\[o{p)dip , = — sm( VNo{p)d(p , (119) 
9 9 

and all other fields vanish. Although the string solutions are gauge equivalent for different 

values of C, the parameter does take on physical meaning when considering multi-string 

configurations in which the value of C is different for different strings Mm . 



Note that the generator (117) can be obtained from by the action of the unbroken 
(electromagnetic) group, 

^ e'CQyig-^Q (120) 

With this in mind, Lepora et al |8^ have classified embedded vortices. The idea is that, 
for a general symmetry breaking G H,th.e Lie algebra of G, G, decomposes naturally into 
a direct sum of the space H of generators of the unbroken subgroup H (the ones associated 
with massless gauge bosons) and the space M of generators associated with massive gauge 
bosons: G = ?{ + M. The action of H on the subspace M further decomposes M into irre- 
ducible subspaces. The classification of embedded vortices is based on this decomposition, 
as we now explain. 

Recall (eq. ^) that finite energy vortices are associated with gauge orbits on the vacuum 
manifold^. Choosing a base point $o in the vacuum manifold, each embedded vortex can 
be associated to a Lie algebra generator which is tangent to the gauge orbit describing 
the asymptotic scalar field configuration of the vortex. The unbroken subgroup H at $o 



"rotates" the various gauge orbits among themselves as in eq. (120). Thus, the action of H 
splits the space of gauge orbits into irreducible subspaces. 

Except for critical values of the coupling constants (which could lead to so-called combina- 
tion vortices), it can be shown |lj, ^ that embedded vortices have to lie entirely in one of 
these irreducible subspaces. If the subspaces have dimension greater than one, then there 
may be a family of gauge equivalent vortices. 

In the GSW model, for instance, the Lie Algebra decomposes into 7i + A^i +A^2 where H is 
spanned by the charge Q, A4i is a one-dimensional subspace spanned by Tz (corresponding 
to the Z-string) and AI2 is a two-dimensional subspace comprising all W-string generators 
Tc- 

Both the W- and the Z-string are embedded string solutions in the GSW model. What 
makes the Z-string more interesting is its unexpected stability properties. It can be shown 
[ [79I that only those vortices lying in one-dimensional subspaces can have a stable semilocal 
limit. Thus, embedded vortices belonging to a family are always unstable. 

Another important difference is that the Z-string is known to terminate on magnetic 
monopoles but this is not true of the W-string. The W-string can terminate without any 
emanating electromagnetic fields since it is entirely within the SU{2) sector of the GSW 
model. 

It is straightforward to embed domain walls in the GSW model. There are no embedded 
monopoles in the GSW model since there is no SU(2) subgroup that is broken to U(l). 
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The gauge orbits are geodesies of a squashed metric on the vacuum manifold 



which is different from the isotropic metric relevant to the scalar sector |81|. 
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6 Electroweak strings in extensions of the GSW model 



Electroweak strings have been discussed in various extensions of tiie GSW model. We de- 
scribe some of tliis work below. We do not, however, discuss extensions in which topological 
strings are produced at the electroweak scale ^ . 



6.1 Two Higgs model 



As discussed in Sec. 5.3, the Z-string is an embedded string in the GSW model. The general 
conditions that enable the embedding are valid even with a more complicated Higgs struc- 
ture. Here we will consider the two Higgs doublet model which is inspired by supersymmetric 
extensions of the GSW model. 

In a two Higgs doublet model, the Higgs structure of the GSW model is doubled so that we 
have scalars $i and $2 and the scalar potential is |Q 

2 2 22 

F($i,$2) = Ai($t<i>i-^)2 + A2($t$2-^)2 + A3[(<i>I<i>i-^) + (<i>^<i>2-^)]' 

+ A4[($I$i)($^$2) - ($l$2)(*^$i)] + A5[Re($l$2) - ^ cose]' 

+ A6[Im($l$2)-^sine]' 

+ A7[Re(<&I$2) - ^ cosC] [Im($I<I>2) - ^ sinC] . (121) 

Here vi and V2 are the respective VEVs of the two doublets, \i are coupling constants and 
the parameter f is a phase 

In polar coordinates, the solution for the two Higgs Z-string is: 

<fi = i.i/i(p)e''^ (^5^ (122) 

*2=i^2/2(p)e''^ (^5) (123) 

(124) 

9z P 

with the profile functions satisfying differential equations similar to the Abelian-Higgs case. 
These have been studied in Ref. |^ where the stability has also been analyzed (also see 

[Ol)- 



6.2 Adjoint Higgs model 



The GSW model with an additional SU{2) field in the adjoint representation, x, is what we 
shall refer to as the "adjoint Higgs model" . The impact of the adjoint field on electroweak 
defects was considered in Ref fzot . 

The bosonic sector of the adjoint Higgs model is: 

L = r„„ + \id^+ ige^W^^U? ~ V{^,x) + Lf (125) 
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where, Te^, is the gradient part of the bosonic sector of the electroweak Lagrangian, i/ is 
the fermionic part of the Lagrangian, efj = eaij {a, i,j — 1, 2, 3) and, 

x) = -A*2*^* - f^lx^ + A2(*^$)^ + Asx* + ax^*^* + 6x • . (126) 

If we impose an additional Z2 symmetry on the Lagrangian under $ —f +$ , x ~^ ~X ■ the 
symmetry is {[SU{2)l x U{1)y]/Z2) x Z2 and we must set 6 = 0. In what follows, we shall 
only consider this case and henceforth ignore the last (cubic) term in the potential. In this 
case, an additional simplification is that the leptons and quarks do not couple to x and so Lf 
is identical to the fermionic Lagrangian of the GSW model. (If the Z2 symmetry is absent, 
the cubic term in the potential is allowed but is constrained to be small by experiment.) 

In a cosmological context, as the universe cools down from high temperatures, if the param- 
eters lie in a certain range |70) there will first be a phase transition in which the adjoint field 
gets a VEV. The VEV of the adjoint will break the SU{2) factor of the high temperature 
symmetry group to U{1). If the VEV of x is along the (0, 0, 1) direction, the generator of 
this U{1) will be and we will denote the unbroken subgroup as U{1)3. So the symmetry 
breaking pattern at this stage is 

([5(7(2) X U{1)y]/Z2) X Z2 ^ ([C/(l)3 X U{1)y]/Z2) x Z2 (127) 

and topological magnetic monopoles will be produced with pure [^(1)3 flux 

At a lower temperature, the doublet field will also get a VEV with the effect, 

([C/(l)3 X U{1)y]/Z2) xZ2^ C/(l)em • (128) 

where, as usual, the electromagnetic charge operator is 

Q - + ^ (129) 

The electromagnetic component (A) from the monopoles is massless but the orthogonal part 
(Z) of the flux is massive and gets confined to a string. This is the Z-string. In addition, the 
breaking of the Z2 factor gives domain walls. 

In the second stage of symmetry breaking, the Z-string is topological and hence is stable. 
The presence of magnetic monopoles from the earlier symmetry breaking means that the 
Z-strings can break by terminating on monopoles. But, as the monopoles form at a higher 
energy scale, their mass is much larger than the energy scale at which strings form and 
which sets the scale for the tension in the string. So the string can only break by instanton 
processes. 

At a yet lower temperature, the VEV of the adjoint turns off. This makes no difference to 
the symmetry structure of the model (apart from restoring the Z2 symmetry and eliminat- 
ing the domain walls) and hence no significant difference to the monopoles connected by 
strings. However, it does affect the stability of the strings since the monopoles are no longer 
topological. 



7 Stability of electroweak strings 
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7.1 Heuristic stability analysis 



As described in |117|, the Z-string goes over into the semilocal string in the Hmit 0^, 7r/2 
and hence the stabihty of the Z-string should match on continuously to that of the semilocal 
string. Therefore we expect that Z-strings should be stable if 9^ is close to 7r/2 and mn < 
mz- 

The stability analysis to certain subsets of perturbations can be carried out much more 
easily than to the completely general perturbations. The subset includes perturbations in 



the Higgs field and W-fields separately. Such analyses may be found in | |117| , [118| , |14| and 
(i) Higgs field perturbations: Perturbations in the Higgs field alone have maximum destabi- 



lizing effect for 9.^, — 7r/4 |14 and, in this case, it is easy to see that the Z-string is unstable. 



Consider the one parameter family of field configurations 

= cosC $o(cosC x) +sm( (130) 

^ji^jO —cos^ Z(o)j(cos^ x) (131) 
where, the string solution is denoted by the subscript, ^ e [0,7r/2] and 



For ^ = the field configuration is the unperturbed Z-string while for ^ = 7r/2 it is the 
vacuum. The energy per unit length of this field configuration can be evaluated and is found 
to be: 

E{£,) ^ cos^ i E{£, ^ Q) . (133) 

Hence the energy per unit length of the string is a monotonically decreasing function of ^ 
and so the string is unstable to decay into the vacuum. 

(ii) Incontractible two spheres: James [^5|, and, Klinkhamer and Olesen fz^ have constructed 
the Z- and W-string solutions by considering incontractible two spheres in the space of elec- 
troweak field configurations in two spatial dimensions. The idea was introduced by Taubes 



[113 1 and was used by Manton to construct the sphaleron [g9|, |73|. The procedure (known 
as the "minimax" procedure) is to construct a set of field configurations that are labelled 
by some parameters /i^. If this set is incontractible in the space of field configurations, then 
there exist (subject to certain assumptions |8^) values of the parameters for which the field 
configuration extremizes the energy functional. For example, Klinkhamer and Olesen Jzi] 
give the following construction for the Z-string in terms of a two parameter (/x, v) family of 
field configurations 

vr/2 < [iJLv] <7r: W = Q , F^O 

$^(l-{l-%)}sin[H)^(j) (134) 

< < 7r/2 : W = -.f{p)G''T- , Y = f{p) sin^ 9.^F^ 

^^h{p)^nu(^^] (135) 



V2 VI 

where, W and Y are Lie algebra valued 1-forms {e.g. W — W^T^da;^), [fif] = max(|/i|, 

F^T" = 2iU-^dU , (136) 
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G°T" = nU[F^T^ + F'^T^ + cos^ e^F^T^]U-^n-^ , (137) 

U (fj,, v,ip) = —i sin ^Ti — i cos /i sin i'T2 — i cos cos v sin 95T3 + cos fi cos cos ipl , (138) 

= C/(^,j/,(^ = 0)"\ (139) 
and the functions f{p) and h{p) satisfy the boundary conditions 

/(O) = h{0) , /(oo) = 1 = h{oo) . (140) 

This set of field configurations labelled by the parameters p^v [— tt, tt] defines an incon- 
tractible two sphere in the space of field configurations. This is seen by considering the 
fields as if they were defined on the three-sphere on which the coordinates are <p, p and v 
and then showing that the field configurations define a topologically non-trivial mapping 
from this to the vacuum manifold which is also an S'^ . Then the minimax procedure says 
that there is an extremum of the energy at some value of the parameters. By inserting the 
field configurations into the energy functional, it can be checked that the extremum occurs 
aX p = Q = V, when the configuration coincides with that of the Z-string. Furthermore, for 
dw < 7''/4, the extremum is a maximum and hence the Z-string is unstable. 

A very similar analysis has been done |6^, Q for the W-string confirming the result [Q 
that it is always unstable. 

(iii) W- condensation: There is also a well-known |l^ instability to perturbations in the 
W-fields alone called "W-condensation" . Application of this instability to the Z-string may 



be found in M, 118, 119, Q]. A heuristic argument goes as follows. 



The energy of a mass to, charge e and spin s particle in a uniform magnetic field B along 
the z-axis is given by: 

E"^ = pI + + {2n + \)eB - 2eB ■ s (141) 

where n — 0,1,2,... labels the Landau levels and Pz is the momentum along the z-axis. 
Now, if s = 1, the right-hand side can be negative for pz — 0, n — provided 

B> — . (142) 

e 

This signals an instability towards the spontaneous creation of spin one particles in sufh- 
ciently strong magnetic fields |Tl| . 

In our case, the magnetic field is a Z-magnetic field and this couples to the spin one W- 
particles. If the string thickness is larger than the Compton wavelength of the W-particles, 
the Z-magnetic field may be considered uniform. Also, the relevant charge in this case is the 
Z-charge of the W-bosons and is gz cos^ 6^ . The constraint that the string be thick so that 
the Z-magnetic field appears uniform and that the charge not be too small means that 6^ 
should be small. Hence the instability towards W-condensation applies for small 0.^. This 
analysis can be performed more quantitatively with the result that there is a relatively 
hard bound sin^ 6w > 0.8 for the string to be stable to W-condensation. 



7.2 Detailed stability analysis 



To analyse the stability of electroweak strings, we perturb the string solution, extract the 
quadratic dependence of the energy on the perturbations and then determine if the energy 
can be lowered by the perturbations by solving a Schrodinger equation. The analysis is quite 
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tedious |6^, ^ |8^, ^ and here we will only outline the main steps. We use the vector 
notation in this section for simplicity. 

The general perturbations of the Z-string are 

{cp^,c^\\,5Z,W%A) (143) 
where, a = 1,2, (j)±^ and 4)\\ are scalar field fluctuations defined by 

\(pNO + h) 

6Z is defined by 

Z = Zno + SZ . (145) 

(The subscript NO means that the field is the unperturbed Nielsen-Olesen solution for the 
string as described in Sec. ^.) The fields W°',A are perturbations since the unperturbed 
values of these fields vanish in the Z-string. 

The perturbations can depend on the z— coordinate and the z—components of the vector 
fields can also be non-zero. However, since the vortex solution has translational invariance 
along the z— direction, it is easy to see that it is sufficient to consider z independent perturba- 
tions and to ignore the z—components of the gauge fields. This follows from the expression 
for the energy resulting from the Lagrangian in eq. (Q) where the relevant z— dependent 
terms in the integrand are: 

lot^G^, + VB.3i^B.3 + (D3<I>)t(i?3<i>) (146) 

and explicitly provide a positive contribution to the energy. Hence we drop all reference to 
the z— coordinate with the understanding that the energy is actually the energy per unit 
length of the string. 

Now we calculate the energy of the perturbed configuration, discarding terms of cubic and 
higher order in the infinitesimal perturbations. We find, 

E=iENo + SENo) + E^+Ec + Ew (147) 

where, E^o is the energy of the Nielsen-Olesen string and SE^o is the energy variation due 
to the perturbations </)|| and SZ. The term E± is due to the perturbation in the upper 
component of the Higgs field, Ec is the cross-term between perturbations in the Higgs and 
gauge fields, while Ew is the contribution from perturbing the gauge fields alone: 



E^^ J cPx + \rf{p - 1)|(/,^|2] , (148) 

Ec = cos6i„ j dPx [$^r'*dj $ - (dj$)^r^$] Wf , (149) 



with Aj defined in [86 



Ew = J iP'x 



ryW^ X . V X i + i|V X VK^ + 7M>2 X Z\ 



where 7 = 5 cos 9^ , 



(150) 



dj = dj - cos{29^)Zj (151) 
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and, the / and Z fields in the above equations are the unperturbed fields of the string. 



The two instabilities discussed in the previous subsection can also be seen in eq. (147). First 
consider perturbations in the Higgs field alone. Then only is relevant. For 0^ = 7r/4, 
dj — dj , and E±_ is the energy of a particle described by the wavefunction (f>± in a purely 
negative potential in two dimensions since < 1 everywhere. It is known that a purely 
negative potential in two dimensions always has a bound state. ^ Hence the energy can be 
lowered by at least one perturbation mode and so the string is unstable when 0m = 7r/4. The 
instability towards W-condensation can be seen in Ew- The term with 7 can be negative 
and its strength is largest for small 0^. Hence W-condensation is most relevant for small 



Returning to the full stability analysis, we first note that the perturbations of the fields that 
make up the string do not couple to the other available perturbations i. e. the perturbations 
in the fields / and v only occur inside the variation SEjsfo- Now, since we know that the 
Nielsen-Olesen string with unit winding number is stable to perturbations for any values 
of the parameters then necessarily, SEnq ^ ^nd the perturbations (/)|| and SZ cannot 
destabilize the vortex. Then, we are justified in ignoring these perturbations and setting 



6Eno = 0. Also we note that the A boson only appears in the last term of eq. ( 150 ) and is 
manifestly positive. So we can set A to zero. 

The remaining perturbations can be expanded in modes: 

- X(p)e™'^ (152) 
for the m*'' mode where m is any integer. For the gauge fields we have, 



{flip) cos{n(p) + flip) sin(n<^)} + ^ {-/ii(p) sininip) + hiip) cosinip)} 
{-hip) sin(n(p) + /2(p) cos(n(^)} ep + ^ {^2(p) cos(nv?) + /i2(p) sin(n(/?)} 



(153) 



(154) 



for the n mode where n is a non-negative integer. 



The most unstable mode is the one with m — and n = 1. This is because these have the 
lowest gradient energy and are the only perturbations that can be non-vanishing ai p — 0. 
Further analysis shows that the string is most unstable to the hi + /12 mode. Hence, we can 
ignore fi, hi — /12 and the barred variables. A considerable amount of algebra then yields: 

SE[x,^+] =27r J dp pix, MO (155) 

where, O is a 2 x 2 matrix differential operator and 

hT + h2 , , 

e+ = . (156) 



Before proceeding further, note that a gauge transformation on the fields does not change 
the energy. However, we have not fixed the gauge in the preceding analysis and hence it 



^''For some potentials though, the wavefunction of the bound state may have 
singular (though integrable) behaviour at the origin and such bound states would 
be inadmissible for us since we require that the perturbations be small. This turns 



out not to be the case for the potential in eq. (148). 
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is possible that some of the remaining perturbations, (x,C+)i might correspond to gauge 
degrees of freedom and may not affect the energy. So we now identify the combination of 
perturbations x ^-nd ^-|_ that arc pure gauge transformations of the string configuration. 

The SU{2) gauge transformation, exp{igip), of an electroweak field configuration leads to 
first order changes in the fields of the form 

(5$ = ig^P'S'o , 6W, = -iuf^ip , (157) 

where Wi — W.°'T°', ip — ifj°'T°-, and the index denotes the unperturbed field and covariant 
derivative. In our analysis above, we have fixed the form of the unperturbed string and so 
we should restrict ourselves to only those gauge transformations that leave the Z-string 
configuration unchanged. (For example, 6^ should only contain an upper component and 
no lower component.) This constrains ■0 to take the form 

^ = sip)( ''7) (158) 



where s{p) is any smooth function. This means that perturbations given by 

are pure gauge perturbations that do not affect the string configuration. Therefore, such 
perturbations cannot contribute to the energy variation and must be annihilated by O. 
Then, in the two-dimensional space of (x, ^+) perturbations, we can choose a basis in which 



one direction is pure gauge and is given by (159) and the other orthogonal direction is the 



direction of physical perturbations. The physical mode can now be written as, 

C(p) = {1- 2 cos' 9^v{p))x{p) + ^^?+(p). (160) 
So now the energy functional reduces to one depending only on C{p)' 

6E[C] =27r f dp p COC (161) 



where O is the differential operator 

p dp \ P+ dp 



= --4-(^4:]+Uip) (162) 



and 



where 



P+ = (1 - 2 cos^ 0^v)' + ^ P f (164) 



g^Tj'f icos^Oy^v'^ , d 2n ^'(l-Scos^ei^i;) 



S+{p) = — — + Pt' 

4 P+{p) dp 



P+ip) 



(165) 



The question of Z-string stability reduces to asking if there are negative eigenvalues lu of 
the Schrodinger equation, 

OC = t^C • (166) 
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Figure 11: The Z-string is stable in the triangular shaded region of parameter space. At 
sin^ 9w = 0.5, the string has a scaling instability. The experimentally allowed parameters 
are also shown. 
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The eigenfunction C(p) must also satisfy the boundary conditions C(/0 = 0) = 1 and C'(0) = 
where prime denotes differentiation with respect to p. In this way the stability analysis 
reduces to a single Schrodinger equation which can be solved numerically. 

The results of the stability analysis are shown in Fig. (jl^) as a plot in parameter space 
(mij/m^, sin^ ^tu), demarcating regions where the Z-string is unstable (that is, where neg- 
ative oj exist) and stable (negative oj do not exist). It is evident that the experimentally 
unconstrained values: sit? 9yj = 0.23 and uiH/rnz > 0.9 lie entirely inside the unstable 
sector. Hence the Z-string in the GSW model is unstable. 

The stability analysis of the Z-string described above leaves open the possibility that the 
string might be stable in some special circumstances such as, the presence of extra scalar 
fields, or a magnetic field background, or fermions. We now describe some circumstances in 
which the Z-string stability has been analyzed. 



7.3 Z-string stability continued 



The stability of Z-strings has been studied in various other circumstances: 

(i) Thermal effects: In pO| the authors examined thermal effects on Z-string stability using 
the high temperature effective potential and found slight modifications to the stability. The 
conclusion is that Z-strings in the GSW model are unstable at high temperatures as well. 
In the same paper, left-right symmetric models were studied and it was found that these 
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could contain stable strings that are similar to the Z-string. 



(ii) Extra scalar fields: It is natural to wonder if the presence of extra scalar fields in 
the model can help provide stability. In p6| the stability was examined in the physically 



motivated two Higgs doublet model with little advantage. In [125| it was shown that an 
extra (globally) charged scalar field could enhance stability. The extra complex scalar field, 
■0, is coupled to the electroweak Higgs by a term |7/'p(f>^'<I> and hence the charges have lower 
energy on the string where <I>^$ '--^ than outside the string where has a non-zero VEV. 
the background of a string and hence a sufficient amount of charge c an s tabilize the string. 
This is exactly as in the case of non-topological solitons or Q-balls |l07| , ^ However, 
scalar global charges attract and this can cause an instability of the charge distribution 



along the string M, 125]. For realistic parameters, stable Z-strings do not seem likely even 



in the presence of extra scalar fields. 

(iii) Adjoint scalar field: A possible variant of the above scheme is that an SU (2) adjoint 



can be included in the GSW model as described in Sec. (6.2). Now, since the Z-string is 
topological within the second symmetry breaking stage in eq. (12S), it is stable. However, 
to be consistent with current experimental data the VEV of the SU{2) adjoint must vanish 
at a lower energy scale. At this stage the Z-string becomes unstable. Hence, in this scheme, 
there could be an epoch in the early universe where Z-strings would be stable. 

(iii) External magnetic field: An interesting possibility was studied by Garriga and Montes 
p6| when they considered the stability of the Z-string placed in an external electromagnetic 
magnetic field of field strength B parallel to the string. First, note that B should be less 
than Be — m^/e, otherwise the vacuum outside the string is unstable to W-condcnsation 
[ pT[ . Then they found that the Z-string could be stable if _B > ^/]3Bc, where P — rrij^/rn^^ 
should be less than 1 for stability of the ambient vacuum. The region of stability for a few 
values of the magnetic field (given hy K = gzB/2m\) is sketched in Fig. |l^. For a certain 
range oi K ^ 0.85, stable Z-strings in the GSW model are still just possible. 

A way to understand the enhanced stability of the Z-string in a magnetic field is to realize 
that the W-condensation instability is due to the interaction of — snvO^A^ + cos6'tuZp 
and . The Z-string itself has a Z magnetic flux. Then the external electromagnetic flux 
can serve to lower the net flux. This reduces the efficiency of W-condensation and 
makes the string more stable. Another viewpoint can be arrived at if we picture the Z- 
string instability to be one in which the string breaks due to the production of a monopole- 
antimonopole pair on the string. If the external magnetic field is oriented in a direction that 
prevents the nucleated magnetic monopoles from accelerating away from each other, it will 
suppress the monopole pair production process, leading to a stabilization of the string for 
sufficiently strong magnetic fields. 

(iv) Fermions: The effect of fermions on the stability of the Z-string has been considered 
in Refs. p7|, |9^ , [76| , |8^ . Naculich |92j found that fermions actually make the Z-string un- 
stable. In |83|] it was argued that this effect of fermions is quite general and also applies 
to situations where the strings form at a low energy scale due to topological reasons but 
can terminate on very massive monopoles formed at a very high energy scale. This most 
likely indicates that the Z-string solution itself should be different from the Nielsen-Olesen 
solution when fermions are included. We shall describe these results in greater detail in Sec. 
^ after discussing fermion zero modes on strings. 

Z-strings have also been considered in the presence of a cold bath of fermions The 
effect of the fermions is to induce an effective Chern-Simons term in the action which then 
leads to a long range magnetic field around the string. 
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Figure 12: The triangular regions depict the parameter range for which the electroweak 
vacuum and the Z-string are both stable in the presence of a uniform external magnetic 
field whose strength is proportional to K. For a range of magnetic field (K ^ 0.85), stable 
strings are possible even with the experimentally constrained parameter values. 
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H Semiclassical stability 



Preskill and Vilenkin |103| have calculated the decay rate of electroweak strings m the region 
of parameter space where they are classically stable. The instability is due to quantum 
tunneling and is calculated by finding the semiclassical rate of nucleation of monopole- 
antimonopole pairs on electroweak strings. The bounce action is found to be 

(167) 

9 as 

where, the strings are classically stable if the ratio of parameters Ooo/os is larger than 1. 
{cLoo/as is the ratio of energy in the magnetic flux when it is spread over an infinite area 
to that if it is confined within the string.) The semiclassical decay probability of the string 
per unit length per unit time is proportional to exp[— 5]. 

The decay rate gets suppressed as we approach the semilocal string {g 0) thus the 
semilocal string is also stable semiclassically. 



8 Superconductivity of electroweak strings 



8.1 Fermion zero modes on the Z-string 



Here we shall consider the fermionic sector of the GSW model in the fixed background of 
the unit winding Z-string for which the solution is given in eq. (82). The Dirac equations 



for a single family of leptons and quarks are obtained from the Lagrangian in Sec. 1.1.2 
These have been solved in the background of a straight Z string in |4^, ^ . The analysis 
is similar to that for U{1) strings since the Z-string is an embedded U{1) string in the 
GSW model (see Sec. 5.3). A discussion of the fermion zero modes in connection with index 
theorems can be found in [iTSl [731 



In polar coordinates with the Z-string along the z-axis, a convenient representation for the 
7 matrices is: 



, (168) 
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(169) 



(Note that the derivative 7^9^^ is given by j*dt+Ydp + l'''d^/p + Y-dz-) Then the electron 
has a zero mode solution 




V^i(p) , en 



ii>i{p) 



(170) 



where, 



P V2 



(171) 
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{q~ l)v u ^ f I 

P V2 



(172) 



In these equations q is the eigenvalue of the operator q defined in eq. (|87|) and denotes the Z- 
charge of the various left-handed fermions. (For the electron, q = cos{29^).) The boundary 
conditions are that tpi and 7/14 should vanish asymptotically. This mean s tha t there is only 
one arbitrary constant of integration in the solution to eqns. ( 171 ) and ( |l72| ). This may be 
taken to be a normalization of ipi and ip/^. 



For the d quark, the solution is the same as in eqns. ( |17C| ), (171) and ( |172| ) except that 
q = 1 — (2/3) sin^ 9yj. For the u quark the solution is: 



UL 




V'2(P) 



iil^aip) 



(173) 



where, 



W2 V2 = -G„— =/V'3 

p \n 

{q + \)v r< ^ f I 

^3 H '03 = -G„— =/V'2 

P V2 



(174) 
(175) 



where, g = -1 + (4/3) sin^ 9^. Note that ( |l7lD , ( |172| ) are related to ( |174| ), ( |175D by q 



The right-hand sides of the neutrino Dirac equations (corresponding to eqns. ( |171| ) and 
( |172| )) vanish since the neutrino is massless. The solutions can be found explicitly in terms 
of the string profile equations in the case when the Higgs boson mass {mu = \/2\rj) equals 



the Z boson mass {mz = gzri/2) 
case are EJ: 



Recall that the string equations in the mn 



f 



yielding the useful relation: 



/ dp— = In 
J P 



mzp 
f 



mz 



(176) 



(177) 



(178) 



where we have included a factor of mz to make the argument of the logarithm dimensionless. 
Now the zero mode profile functions for the massless fermions are: 



; 3/ 



3/2 f mzp\ 
^ \ f J 



tpi = C4TO2 



3/2/ rnzp\ 
^ \ f J 



9-1 



(179) 



where, Ci and C4 are independent constants that can be chosen to normalize the left- and 



right-handed fermion states and the spinors are given in (170). The boundary condition 
that the left-handed fermion wavcfunction should vanish at infinity is only satisfied if q > 0. 



Hence (179) can only give a valid solution for q > for the left-handed fermion. If we also 
require normalizability, we need q > 1. (Note that there is no singularity at p = because 
/ oc p when p ~ 0.) If we have a left-handed fermion with q < —1, the correct equations 
to use are the equations corresponding to the up quark equations given in (174) and (175) 
and these are solved by letting q —q in (179). In this case, the spinors are given in (173). 



For the electroweak neutrino, the right-handed component is absent and q = — 1. This means 
that the neutrino has the same spinor structure as the left-handed up quark and the solution 
is that in (179) with q replaced by -1-1. Therefore the wave function falls off as 1/p and 
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Figure 13: The direction of propagation of quark and lepton zero modes on the Z-string. 
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the state is strictly not normahzable - the normahzation integral diverges logarithmicaUy. 
However, depending on the physical situation, one could be justified in imposing a cut-off. 
For example, when considering closed loops of string, the cutoff is given by the radius of 
the loop. 



Next we give the explicit solutions to the Dirac equations in (171) and (172) in the case 
when the fermion mass {nif = hrj/y/2) is equal to the scalar mass which is also equal to 
the vector mass. This so-called "super-Bogomolnyi" limit is not realized in the GSW model 
but may be of interest in other situations (for example, in supersymmetric models). Then, 
if the charge on the left-handed fermion vanishes ((7 = 0), the solution can be verified to be: 

V'4(p) = 2iVmf :M(l-^;(p)) (181) 
P 

where A'^ is a dimensionless normalization factor. For the same set of parameters, the solution 
for the up quark equations can be written by using the transformation q ^ —q in the 
above solutions. Further, this solution can also be derived using supersymmetry arguments 



[126, \m. 



The left-handed fermion wave-functions found above can be multiplied by a phase factor 
eyip[i{Ept — pz) and the resulting wave-function will still solve the Dirac equations provided 

Ep = e,p (182) 

where, i labels the fermions, and, 

ei/ = +1 = e„ , Ce = -1 = Ed • (183) 



In other words, vl and u travel parallel to the string flux while e and d travel anti-parallel 
to the string flux. 

We should mention that the picture of quarks travelling along the Z-string may be inaccurate 
since QCD effects have been totally ignored. At the present time it is not known if the strong 
forces of QCD will confine the quarks on the string into mesons and baryons (for example, 
pions and protons). Further, the electromagnetic interactions of the particles on the string 
might lead to bound states of electrons and protons on the string. This would imply a 
picture where hydrogen (and other) atoms are the fundamental entities that live on the 
string. 
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8.2 Stability of Z- string with fermion zero modes 



In Fig. |lj we show the effect that perturbations of order e in the Z-string fields have on the 
fermion (u and d quarks) zero modes. The zero momentum modes acquire an 0(e) mass 
while the non-zero momentum modes get an O(e^) mass. For the perturbation analysis to 
make sense, we require that the u and d quark zero momentum modes are either both filled 
or both empty. In that case, the 0(e) terms in the variation in the energy will cancel and 
we will be left with something that is O(e^). In fact, 

,2 N T 

Ai? = --KpL^- (184) 
fe=i 

where mi is a matrix element having to do with the interactions of the u and d quarks, 
L — !■ cx) is the length of the string on which periodic boundary conditions have been imposed, 
and — > oo is a cut-off on the energy levels which are labeled by k. The crucial piece of 
this formula is the minus sign which shows that the energy of the string is lowered due to 
perturbations fo^ . 

In Ref. it was argued that an identical calculation could be done for any classically 
stable string that could terminate on (supermassive) magnetic monopoles. However, in the 
low energy theory, the strings are effectively topological and hence, it seems unlikely that 
fermions can lead to an instability. This suggests that the bosonic string configuration gets 
modified by the fermions and the stability analysis around the Nielsen-Olesen solution may 
be inappropriate. 

So far, the stability analysis with fermions presented here only considered the zero modes 
and ignored the infinitely many massive fermion modes. Very recently, Groves and Perkins 
have analysed the full spectrum of massless and massive fermionic modes in the back- 
ground of the electroweak string. They then calculate the effect of the Dirac sea on the 
stability of electroweak strings by calculating the renormalised energy shift of the Dirac 
sea when a Z-string is perturbed by introducing a non-zero upper component to the Higgs 
doublet. This energy shift is negative and so destabilises the string, but it is small, leading 
them to conclude that if positive energy fermionic states are populated, it is conceivable 
that the total fermionic contribution could be to stabilise the string. This work is still in 
progress. In the meantime, the stability of Z-strings remains an open question. 



8.3 Scattering of fermions off electroweak strings 



The elastic scattering of fermions off semilocal and electroweak strings has been considered 

in n, gg. 

The main feature of the cross section is that the scattering violates helicity |Q. It is 
straightforward to show that the helicity operator S • 11, where = e'^'=^*--y-' is the spin 
operator and 11* are the canonical momenta, does not commute with the hamiltonian. If 
i"-^ = (^+, (fP), the commutator is proportional to {D(jf) terms. Consider for a moment the 
usual representation of Dirac matrices, 

^-0 ;) ^-(°.-o') -(o'-°o <-) 
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Figure 14; The effect of perturbations of the Z-string on fermion zero modes. 
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Then, for an incoming electron, one finds 



[H,S.^]=^;^(^ ^,.^^,.^0 ), (186) 

where h is the Yukawa couphng and {DjcjP) is given in eq. (|8^). Therefore hehcity- violating 
processes can take place in the core of the string. 

A preliminary calculation by Ganoulis in ref. showed that, for an incoming plane wave, 
the dominant mode of scattering gives identical cross sections for positive and negative 
helicity scattered states. More precisely, for an incoming electron plane wave of momentum 
k, energy w and positive helicity it was found that, to leading order. 



dk 



1 / ui — ^ 



k \ 2u} 



sm^TrqR) (187) 



k"^ + mg, qh is the Z-charge of the right fermion field, given in eq. ( p7| ) (recall 
that right and left fermion fields have different Z-charges, = gi, ± 1). 

A more detailed calculation was done by Davis, Martin and Ganoulis |Q , and later extended 
by Lo using a 'top hat' model 

^W^V;/V2 r>R ""^'^ = \2/g. r>R ' ^^^^^ 

which is expected to be a reasonable approximation since the scattering cross section in 



the case of cosmic strings has been shown to be insensitive to the core model |101]. Note 
that there is a discontinuous jump in the fermion mass and string flux; however the wave 
functions are matched so that they are continuous at r = i?. Note that the left and right 
fields decouple in the core of the string, so helicity violating processes are concentrated at 

r = R. 

The authors of ^ confirmed that, in the massive case, there are helicity-conserving 
and helicity-flip scattering cross sections of equal magnitude. The latter goes to zero in the 
massless limit (in that case, the left and right fields decouple, and no helicity violation is 
possible), suggesting that helicity violation may be stronger at low energies. For "fractional 
string flux" {i.e. for fractional q) the cross section is of a modified Aharonov-Bohm form, and 
independent of string radius. For integer q it is of Everett form ||40|| (the strong interaction 
cross section is suppressed by a logarithmic term). 

Another interesting feature has to do with the amplification of the fermionic wave function 
in the core of the string. Lo has remarked that there is a regime in which the scattering 
cross section for electroweak strings is much less sensitive to the fermion charge (that is, to 
sin^^u)) than for cosmic strings. In contrast with, e.g., baryon number violating processes, 
which show maximal enhancement only for discrete values of the fractional flux, the helicity 
violating cross section for electroweak strings in the regime fc ~ m, kR << 1 shows a 
plateau for < sin^ < 1/2 where amplification is maximal and the cross section becomes 
of order wj^- This can be traced back to the asymmetry between left and right fields; 
while the wave function amplification is a universal feature, different components of the 
fermionic wave function acquire different amplification factors in such a way that the total 
enhancement of the cross section is approximately independent of the fermionic charge, q 
(or, equivalently, of sin^O^). 

Elastic scattering is independent of the string radius for both electroweak and semilocal 
strings (for integral flux there is only a mild dependence on the radius coming from the 
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logarithmic suppresion factor in the Everett cross section). Since the cross section is hke 
that of U{1) strings, we would expect electroweak and semilocal strings to interact with the 
surrounding plasma in a way that is analogous to topological strings. 



9 Electroweak strings and baryon number 



As first shown by Adler [gj, and, Bell and Jackiw [Q, currents that are conserved in a 
classical field theory may not be conserved on quantization of the theory. In the GSW model, 
one such current is the baryon number current and the anomalous current conservation 
equation is: 



327r2 



(189) 



where jg is the expectation value of the baryon number current operator bs '■ ip^^ip : 
where the sum is over all the species of fermions labeled by s, V' is the fermion spinor and 
hs is the baryon number for species s and the operator product is normal ordered. Also, Np 
denotes the number of families, and tilde the dual of the field strengths. 

The anomaly equation can be integrated over all space leading to 

AQs = ^ y dtSx[~g^W;,W-^- + g'^Y^^Y^-] ^ AQcs- (190) 



with, 



32^2 / « ^^yfe 



(191) 



Here, A(-) denotes the difference of the quantities evaluated at two different times, Qb is 
the baryonic charge and the surface currents and integrals at infinity are assumed to vanish. 
Qcs is called the Chern-Simons, or topological, charge and can be evaluated if we know 



the gauge fields. The left-hand side of eq. (190) evaluates the baryon number by counting 
the fermions directly. We describe the evaluation of both the right- and left-hand side for 
fermions on certain configurations of Z-strings in the following subsections. Finally, in Sec. 
9.4 we briefiy comment on possible applications to cosmology. 



9.1 Chern-Simons or topological charge 



We will be interested in the Chern-Simons charge contained in configurations of Z-strings. 
Then, we set all the gauge fields but for the Z-field to zero in the expression for the Chern- 
Simons charge, yielding 

2 P 

Qcs = Nf:^ cos(20^) / d^xZ ■ Bz (192) 
327r^ J 

where, Bz denotes the magnetic field in the Z gauge field: — ^^•'''djZk- 

The terms on the right-hand side have a simple interpretation in terms of a concept called 
"helicity" in fiuid dynamics [l8| . Essentially, if a fluid flows with velocity v and vorticity 
uj — W X V, then the helicity is defined as: 

d^xv-uj (193) 
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Since the helicity measures the velocity flow along the direction of vorticity, it measures 
the corkscrew motion (or twisting) of the fluid flow. A direct analog is defined for magnetic 
fields: 

(fxA-B (194) 



which is of the same form as the terms appearing in ( 192 ) . Hence the Chern-Simons charge 
measures the twisting of the magnetic lines of force. The helicity associated with the Z field 
alone is given by: 



(FxZ ■ Bz 



(195) 



If we think in terms of flux tubes of Z magnetic field, Hz measures the sum of the link and 
twist number of these tubes: 

Hz = Lz+Tz . (196) 

For a pair of unit winding Z flux tubes that are linked once as shown in Fig. |l^ the helicity 
is: 



2Fi 



(197) 



where, Fz is the magnetic flux in each of the two tubes Note that the helicity is positive 
for the strings shown in Fig. |l^. If we reversed the direction of the flux in one of the loops, 
the magnitude of Hz would be the same but the sign would change. For the Z-string, we 
also know that 

47r 

Fz = — (198) 

9z 



and so eq. (192) yields |123|: 



Qcs^NFC0s{2e^) 



(199) 



9.2 Baryonic charge in fermions 



The baryon number associated with linked loops of Z-string has been evaluated in Ref. 
by studying the fermionic zero modes on such loops. This corresponds to evaluating the 



left-hand side of eq. (190) directly in terms of the fermions that carry baryon number. The 
calculation involves adding the baryonic charges of the infinite Dirac sea of fermions living 
on the string together with zeta function regularization. 

To understand why the linking of loops leads to non-trivial effects, note that the quarks and 
leptons have a non-trivial Aharanov-Bohm interaction with the Z-string. So the Dirac sea 
of fermions on a loop in Fig. ^ is affected by the Z-ffux in the second loop. This shifts the 
level of the Dirac sea in the ground state leading to non-trivial baryonic and other charges. 

Instead of considering the linked loops as shown in Fig. |l5| it is simpler to consider a large 
circular loop of radius a oo in the a;?;-plane threaded by n straight infinite strings along 
the z-axis (Fig. [l6| ). Then the fermionic wave- functions take the form: 

= e-'(^''*-P'^VL^(?-) , i^fl = e-*(^''*-(P-"/")")i/^^"V) (200) 
where the functions with superscript (0) are the zero mode profile functions described in Sec. 



.1 



and CT is a coordinate along the length of the circular loop. From these wavefunctions, 
the dispersion relation for a zero mode fermion on the circular loop is 

tuk = e,{k - qZ) . (201) 
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Figure 16: A circular Z-string loop of radius a threaded by n Z-strings. 



n Z-strings 




where, q is the Z-charge of the fermion, ei is defined in cq. (182), uj is related to the energy 
E hy uj = aE, and k to the momentum p by k = ap € E. Z is the component of the gauge 
field along the circular loop multiplied by a and is given by, 

Z= — . (202) 

9z 

The crucial property of the dispersion relation is that, if there is an Aharanov-Bohm inter- 
action between the Z-string and the fermion, ujk cannot be zero for any value of k since k 
is an integer but qZ is not. 

The Z- , A- and baryon number (B) charges of the leptons and quarks are shown in Table 
|l|. Note that we use 2qz/gz to denote the Z-charge and this is identical to the eigenvalue of 
the operator q defined in eq. |8^ and also to q used in the previous section. 

The energy of the fermions is found by summing over the negative frequencies - that is, the 
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Table 1: Summary of , electric and baryonic charges for the leptons and quarks. The 
charges qz are for the left-handed fermions and s" = sin" 9^ ■ 







e 


d 


u 




-1 


1 - 2s" 


1 - 2s"/3 


-1 + 4s"/3 







-1 


-1/3 


2/3 










1/3 


1/3 



Dirac sea - and so the energy E due to a single fermion species is: 

n ^ — ^ n ^ — ^ 



(203) 



k—kf 



where, kp denotes the Fermi level - the value of k for the highest filled state. Therefore we 
need to sum a series of the type: 



J2ik-qZ) = Xl(fc + kF- qZ) 



k—kp k—0 

The sum is found using zeta function regularization 

1 1 



S ^a-'^.kp-qZ) 



[kp — qZ){kp — qZ — 1) 



(204) 



(205) 



12 2 

With this result, the energy contribution from the z*'' species of fermions takes the form: 



E,. 



1 1 

24a ^ 2a 



Si) 



q^Z+^ 



1 1 , 

1 

24a 2a ' 



(206) 



Adding the contributions due to different members of a single fermion family, we get 



2a£: = Kl + Kl + 3X? + iKl 



(207) 



Next we can calculate the angular momentum of the fermions in the circular loop back- 
ground. The system has rotational symmetry about the z— axis and this enables us to 
define the generalized angular momentum operator as the operator that annihilates the 
background field configuration |69] : 



where. 



dip 



(208) 



(209) 



Sz is the spin operator, and, the isospin operator is given in terms of the U (1) (hypercharge) 
and SU (2) charges - qi and q2 respectively - of the field in question: 



292 \y3 



2cn 
9' 



(210) 



The isopin operator acts via a commutator bracket on the gauge fields and by ordinary 
matrix multiplication on the Higgs field and fermion doublets. 



We are interested in the angular momentum of the chiral fermions on the circular loop which 
lies entirely in the xy-plane. The fermions in the zero modes therefore have Sz — 0. (The 
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spin of the fcrmions is oriented along their momenta which hes in the xy-plane.) The action 
of L2 is found by acting on the fermion wave-functions such as in eq. ( ^OOD (remembering to 
let n — + ~n for the neutrino and up quark). The action of Iz is found by using the charges 
of the fcrmions given in the GSW model defined in Sec. [TT^ . We then find: 



ul 
dL 



(211) 



M^CR = k^'^en , M,UR = + ji/i)uR , M,dR = (fc^"^) - 2n/i)dR (212) 
where the fc'^*-' are defined above eqn. (|20l|). Now summing over states, as in the case of 



the energy, we find the total generalized angular momentum of the fermions on the circular 
loop: 



1 




2 


1 


\(-) 


1 


2 


3 


2 






2 




2 




2 



id) 2n 1 



1 2 



n 
3 



1 2 



(213) 



Note that though the gauge fields do not enter explicitly in the generalized angular momen- 
tum, they do play a role in determining the angular momentum of the ground state through 
the values of the Fermi levels. 

The calculation of the electromagnetic and baryonic charges and currents on the linked 
loops is similar but has a subtlety. To find the total charge, a sum over the charges in all 
filled states must be done. This leads to a series of the kind: 



To regularize the divergence of the series, it is written as 

{k - qZ f . 



Sn = lim V "'^^^ 



(214) 



(215) 



k—kji 



The subtlety is that the gauge invariant combination k — qZ is used as a summand rather 
than k or some other gauge non-invariant expression pO| . Once again zeta function regu- 
larization is used to get: 



Sg = ^(/c + kF- qZf = C(0, kp -qZ) = - 



k=0 



With this result, the contribution to the charge due to fermion i is: 



kf-q.Z+'-^ 



(216) 



(217) 



where, qi is the charge carried by the i*-^ fermion of the kind that we wish to calculate. 
(Note that qi can represent any charge - electric, baryonic etc. - and is, in general, different 
from the Z-charge qi.) 



The currents along the string are given by ip^^tp where 7^ is given in eq. ( [1691) . This gives 



Ji — 



2'Ka 



(218) 



By adding the contributions due to each variety of fermion, expressions for the energy, 
angular momentum, charges and currents for one loop threaded by n have been found in 
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Figure 17: The energy of the ground state of linked loops versus x = 2nsin^ 0^/3. 
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Table 2: Expressions for the energy, generalized angular momentum, charges and currents 

in terms of a; = 2nsiv? 9w/<!>- We have omitted the multiplicative factor Np in all the 
expressions for convenience. 

a: G (0,1/3) (1/3,1/2) (1/2,2/3) (27M) 

12a;^ -Qx + l 12a;^ - 9a; + 2 12a;^ - 15a; + 5 12a;^ - 18a; + 7 

M n-l 2-n 

Qa/c 0-1+1 

B -3a; + 1 -3x + 1 -3x + 2 -3a; + 2 

2TTaJA/e -8a; + 2 -8a; + 3 -8x + 5 -8a; + 6 

2'KaJB —X —X 1 — X 1 — X 



[ [47| . These results are reproduced in Table It is reassuring to note that in the ground 
state, the baryon number of the single loop is given by nNp cos 20w in agreement with the 
calculation of the Chern-Simons number. 

The energy of the fermionic ground state shows a complicated dependence on x as is demon- 
strated in Fig. |l^. Note that E{x) does not have a monotonic dependence on x and the 
energy of strings that are linked n times bears no simple relation to those linked m times. 
In particular, the energy does not continue to decrease as we consider strings that have 
higher linkage. The lowest energy possible, however, is when a; = 1/4 and for n — 1, this 
corresponds to sin^ 9^ = 3/8, which is also the value set by Grand Unified models. It is not 
clear if this is simply a coincidence or if there is some deeper underlying reason p39| . 



9.3 Dumbells 



In his 1977 paper, Nambu discussed the possible occurrence of electroweak monopoles and 
strings in particle accelerators. There are two issues in this discussion: the first is the pro- 
duction crossection of solitonic states in particle collisions, and the second is the signatures 
of such states if they are indeed produced in an accelerator. The answer to the first question 
is not known though it is widely believed that the process is suppressed not only by the 
large amount of energy required but also due to the coherence of the solitonic state. The 
second question was addressed by Nambu |Q and he estimated the energy and lifetime of 
electroweak strings that may be possible to detect in accelerators. 

To find the energy of a Z-string segment, Nambu treated the monopoles at the ends as hollow 
spheres of radius R inside which all fields vanish. A straightforward variational calculation 
in units of 77 « 246 GeV then gives the monopole mass 



M=^sin^/^^../^ (219) 

and radius 

R ^ l^K. (220) 
V mnrnw 

The string segment is approximated by a cylindrical tube with uniform Z magnetic flux 
with all other fields vanishing. This gives 



r^A—\ (221) 
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for the core radius and string tension, respectively. 



Now, if the monopoles are a distance / apart, the total energy of the system is 

E = 2M - — +tI (222) 

which is clearly minimised by Z = i.e. the string can minimize its energy by collapsing. 
The tendency to collapse can be countered by a centrifugal barrier if the string segment 
("dumbell") is rotating fast enough about a perpendicular axis. The energy and angular 
momentum of a relativistic dumbell has been estimated by Nambu to be: 

E^^ttIt L^-ttI^t (223) 
8 

where, 

— = ^ (224) 

with V ^ 1 being the velocity of the poles. The expressions for E and L imply the existence 
of asymptotic Regge trajectories, 

L - a'^E^ (225) 

with slope 



, 1 ( mz 
0^0=0 

which, if found, would be a signature of dumbells 



TeV^^ . (226) 



The orbiting poles at the ends of the rotating dumbell will radiate electromagnetically and 
this energy loss provides an upper bound to the lifetime of the configuration. An estimate of 
the radiated power from the analysis of synchrotron radiation in classical electrodynamics 
(eg. see Q) gives 



_ X 137 (^—j sin^a^ (227) 

Therefore the decay width T = P/E is given by 

E 

r ~ — (228) 
and for large angular momentum, can lead to significant lifetimes (compared to E^^). 

To obtain numerical estimates, note that the above estimates are valid only if the dumbell 
length is much greater than the width of the Z-string. This imposes a lower bound on the 
angular momentum: 

L » I X 137 sin^ 9^, cos^ 9^ ~ 36 (229) 



Using the relation between the energy and the angular momentum, such an object has 
E » <o{mHlmzfl'^ TeV. 

The estimates above assume that the lifetime of the dumbell is dictated by the energy 
emission in photons. In reality, there are other decay channels as well, though it is likely 
that these will be comparitively suppressed since the photon is the only massless boson 
present in the system. The dumbell can also decay by fragmenting due to field-theoretic 
instabilities of the kind discussed in Sec. 0. These may be suppressed due to the finite size 
of the dumbell, and as Nambu points out, due to the angular momentum of the dumbell |^ 
A careful analysis of these factors has not yet been performed and is a vital open problem 
that may become experimentally relevant with the next generation of accelerators. 

^^In the stability analysis for a finite piece of string of length L, the eigenvalues of 
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9-4 Possible cosmological applications 



The role of electroweak strings in cosmology depends on their abundance during and after 
the electroweak phase transition. If this abundance is negligible, electroweak strings may at 



best only be relevant in future accelerator experiments (see Sec. 9.3). If, however, there is 
a cosmological epoch during which segments and loops of electroweak strings were present, 
they could impact on two observational consequences: the first is the presence of a primordial 
magnetic field, and the second is the generation of a cosmological baryon number. What 
is perhaps most remarkable is that the two consequences might be related - the baryonic 
density of the universe would be related to the helicity of the primordial magnetic field 



119| , |106 | 



(i) Primordial magnetic fields: A gas of electroweak segments is necessarily accompanied 
by a gas of electroweak monopoles. The eventual collapse and disappearance of electroweak 
strings removes all the electroweak monopoles but the long range magnetic field emanating 
from the monopoles is expected to remain trapped in the cosmological plasma since that is 
a very good electrical conductor. This will then lead to a residual primordial magnetic field 
in the present universe. 

A quantitative estimate of the resulting primordial magnetic field cannot be made with 
confidence but a dimensional estimate is possible. An estimate for the average flux through 
an area — N"^ /T"^ , where is a dimensionless number that relates the lengt h scale o f 



interest, L, to the cosmological thermal correlation length T ^, was obtained in |116, 119| 



and then translated into the average magnetic field through that area. The result is: 

B\area ~ T^/A . (230) 

(Magneto-hydrodynamical considerations provide a lower bound ^ 10^^ cms on L at the 
present epoch.) It is important to remember that the above is an areal (i.e. flux) average, 
defined by [|§ 

jJdS-B^ ^ (231) 

where the surface integral is over an area A and (•) denotes ensemble averaging. 

(ii) Baryon number: A gas of electroweak string segments and loops would, in general, con- 
tain some helicity density of the Z-field. When the electroweak strings eventually annihilate. 



it is possible that the helicity gets converted into baryon number 123 , 11E |. However, in 
Ref. [^l], ^ it is argued that fractional quantum numbers of a soliton are unrelated to 
the number of particles produced when the soliton decays. Instead, only the change in the 
winding of the Higgs field in a process that starts out in the vacuum and ends up in the 
vacuum can be related to the particle number. This would imply that we would have to 
consider the formation of electroweak strings together with their decay before we can find 
the resulting baryon number. Such a calculation has not yet been attempted. 

An interesting question is to consider what happens to the helicity in the Z-ficld after the 
strings disappear. One possibility is that the helicity gets transferred to a frozen- in residual 
magnetic field after the strings have decayed. To see this, consider a linked pair of loops 



the stability equation are shifted by a contribution of order vr^/L^ with respect to 
the infinitely long case, thus for sufficiently short segments the radial decay mode 
could become stable. Longitudinal collapse might then be stabihzed by rotation, as 
explained above 
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as in Fig. |T^. The strings can break by nucleating monopole-antimonopole pairs, and then 
the string segments can shrink, finally leading to monopole annihilation. If this process 
happens in the early universe, the loops will be surrounded by the ambient plasma which 
will freeze-in the magnetic field lines. Hence, after the strings have disappeared, we will be 
left with a linked pair of magnetic field lines. In other words, the original helicity in the 
Z-field has been transferred to helicity in the A-field. This argument relies on the freezing-in 
of the magnetic field emanating from the monopoles and in the real setting the physics can 
be much more complicated. However, a connection between the baryon abundance of the 
universe and the properties of a primordial magnetic field seems tantalizing. 



Stable strings at the electroweak scale: If in more exotic models, strings at the electroweak 
scale were stable and had the superconducting properties discussed above, they could be 
responsible for baryogenesis [|l3[ and the presence of primary antiprotons in cosmic rays 



112 1 . The production of antiprotons follows on realizing that any strings tangled in the 



galactic plasma would be moving across the galactic magnetic field. In the rest frame of 
the string, the changing magnetic field causes an electric field along the string according to 
Faraday's law. The electric field along the string raises the levels of the u- and d-quark Dirac 
seas (see Fig. |^), as well as the electron Dirac sea (not shown in the figure). This means 
that the electric field produces quarks and leptons on the string. The electric charges of the 
particles are in the ratio e : u : d :: —I : +2/3 : — 1/3 and the rate of production of these 
particles due to the applied electric field is proportional to the charges. Furthermore, the 
quarks come in three colors and so for every electron that is produced, 3x2/3 = 2 u-quarks 
and 3 x 1/3 — 1 d-quark are also produced. As a result, the net electric charge produced is 
1 X (-1) + 2x (2/3) 1 X (-1/3) = 0. However, net baryon number 2 x (1/3) -f 1 x (1/3) = 1 
is produced because the quarks carry baryonic charge 1/3 while the baryonic charges of the 
leptons vanish. Depending on the orientation of the string, either baryons or antibaryons 
will be produced. Some of these would then be emitted from the string and would arrive on 
earth as cosmic rays. 



Formation of strings in the electroweak phase transition. Early attempts to understand the 
formation rates of electroweak strings were made in |119 based on the statistical mechanics 



of strings. The estimates indicate that a density of strings will be formed immediately after 
the phase transition. However, the application of string statistical mechanics to electroweak 
strings may not be justified and so other avenues of investigation are needed. An alternative 
approach to study electroweak string formation was taken by Nagasawa and Yokoyama |0 . 
They assumed a thermal distribution of scalar field values and gradients, and estimated 
the probability of obtaining a string-like scalar field configuration. The conclusion was that 
electroweak vortex formation in a thermal system is totally negligible. One possible caveat is 
that the technique used in ignores the effect of gauge fields, which we know are significant 
in the formation of related objects such as semilocal strings. In |109|, Saffin and Copeland 
have evolved the classical equations of motion to study the formation of electroweak strings, 
and they found the presence of the gauge fields led to larger string densities than one would 
have inferred from the scalar fields alone, at least when sin^ 9^^ — 0. However, this study 
does not directly address the question of string formation in a phase transition because no 
measure has been placed on the choice of initial conditions and their choice may be too 
restrictive. Most recently, a promising development has taken place pGf - calculations in 
lattice gauge theory have been done to study the electroweak phase transition and there is 
evidence that electroweak strings will form. Further studies along these lines will provide 
important and quantitative insight into the formation of electroweak strings. 



Using the results on the formation of semilocal strings, we can gain some intuition about 
the formation of electroweak strings in the region of parameter space close to the semilocal 
limit (the region of stability in Fig. |ll|) . We have seen that semilocal strings with (3 < I have 
a non-zero formation rate, increasing as /3 — ^ 0. Initially short segments of string are seen 
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Figure 18: The dispersion relations for the u and d quark zero modes are shown. The filled 
states are denoted by solid circles while dashes denote unfilled states. For convenience, 
periodic boundary conditions are assumed along the string and so the momentum takes on 
discrete values. 




to grow and join nearby ones because this reduces the gradient energy at the ends of the 
strings. The ends of electroweak strings are proper magnetic monopoles, and therefore the 
scalar gradients are cancelled much more efficiently by the gauge fields, but as sin 
the cores of the monopoles get larger and larger, and they could begin to overlap with 
nearby monopoles, so it is possible that short segments of electroweak string will also grow 
into longer ones. 



10 Electroweak strings and the sphaleron 



The sphaleron is a classical solution in the GSW model that carries baryon number Nf/2, 

|. For 6m = 0, the asymptotic form of 
the sphaleron Higgs field is: 

'^-''-(sin?!-) ■ (232) 

in which $ should be replaced 
by ^sph- (Note that the hypercharge gauge field vanishes for 6*^ = 0.) Inside the sphaleron, 
the Higgs field vanishes at one point. The sphaleron also has a magnetic dipolc moment 
that has been evaluated for small values of . The reason that the sphaleron is important 
for particle physics is that its energy defines the minimum energy required for the classical 
violation of baryon number in the GSW model. 

As has already been described in Sec. ^, non-trivial baryon number can be associated with 
linked and twisted segments of electroweak string. Further, for specific values of the link 
and twist, the baryon number of a configuration of Z-strings can also be Np/^. This raises 
the question: are sphalerons related to Z-string segments? 

An early paper to draw a connection between the various solutions in the GSW model is Ref. 

, however, a direct correspondence between the field configuration 
of the Z-string and the sphaleron was made. 



where Np is the number of fermion families 




while the gauge fields continue to be given by eq. (|lOd)-(101 



In 1123 m 119 
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10.1 Content of the sphaleron 



In [p9[ Hindmarsh and James evaluated the magnetic charge density and current density 
within the sphaleron. A subtlety in this calculation is that there is no unique definition of 
the electromagnetic field when the Higgs field is not everywhere in the vacuum. The choice 
adopted in |59| (and also the choice in this review) is 

F^^J^ = sin 61.^ W^^. n'' + cos 9^ Y,j . (233) 

The evaluation of the magnetic charge density (which is proportional to the divergence of 
the magnetic field strength) clearly shows that the sphaleron contains a region with positive 
magnetic charge density and a region with negative magnetic charge density. Furthermore, 
the total charge in, say, the positive charge region agrees with the magnetic charge of a 
monopole. In addition, there is a flux of Z magnetic field connecting the two hemispheres. 
This would seem to confirm that the sphaleron consists of a Z-string segment. However, 
this is not the full picture. In addition to the string segment, Hindmarsh and James find 
that the electric current is non-zero in the equatorial region and is in the azimuthal {e^p) 
direction. 



10.2 From Z- strings to the sphaleron 



_ , cos(e/2) . 



The scalar field configuration for a finite segment of Z-string was given in Sec. 3.3 

cos(e/2) 
_sin(e/2) e'f 

where, 

cos Q = cos 9,n — cos 9ffi + 1 (235) 

and the angles 0„i and Om are measured from the monopole and antimonopole respectively, 
as shown in Fig. n9. 



It is straightforward to check that (234) yields the monopole field configuration close to 
the monopole {9m ~* 0) and the antimonopole configuration close to the antimonopole 
{9m — '"■)• It also yields a string singularity along the straight line joining the monopole 
and antimonopole {9m — t^, 9m — 0). However, there are other Higgs field configurations 
that also describe monopoles and antimonopoles: 

_ f cos{em/2) \ ff> f sin(0^/2) \ , 

[sm{9m/2) e^^ ) ' ^™ " ^ [cos{9r^/2) e^^ ) ' ^^^^^^ 

Next consider the Higgs field configuration: 

,._f sm{em/2)sm{9m/2)e'^ + cos{9m/2)cos{9m/'2) \ , . 

^ramKl) - l^sin(0„,/2)cos(0^/2)e'^ - cos(0,„/2) sin(0^/2)e'(^-^) ) ^^-^ > 



together with the gauge fields given by eq. ( |100D -(101) with $ replaced by ^mfh{l)- When 
we take the limit 9m — > we find the monopole configuration (with 7 = 0) and when we 
take 9m — > tt the configuration is that of an antimonopole (with arbitrary 7) provided we 
perform the spatial rotation ^ -f 7. Note that the asymptotic gauge fields agree since 



these are determined by the Higgs field. The monopole and antimonopole in (237) also have 



the usual string singularity joining them. This means that the configuration in eq. ( 237 ) 



describes a monopole and antimonopole pair that are joined by a Z-string segment that is 
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Figure 19: Definition of the coordinate angles 9m and Ofn- The azimuthal angle, is not 
shown. 




twisted by an angle 7. The Chern-Simons number of one such segment can be calculated 



123 1 and is 



Qcs ^ Nf cos20^ ^ . (238) 

If 7 = tt/ cos(26'tu) then the Chern-Simons number of the twisted segment of string is Nf/2 
and is precisely that of the sphaleron. 

Given that the segment with twist tt/ cos(20.u,) has Chern-Simons number equal to that of 
the sphaleron, it is natural to ask if some deformation of it will yield the sphaleron. This 
deformation is not hard to guess for the 9^ — case. In this case, if we let the segment size 



shrink to zero, we have 9m — Om — and the Higgs field configuration of eq. (237) gives: 



This is exactly the scalar field configuration of the sphaleron for 9^ — (eq. (^32|)). Note 



that the asymptotic gauge fields continue to be given by eq. (IOC)- (101) and satisfy the 



requirement that the covariant derivatives of the Higgs field vanish. 



Encouraged by this successful connection in the 9^, = Q case, it was conjectured in |122, 119 
that the sphaleron can also be obtained by collapsing a twisted segment of Z-string with 
Chern-Simons number Np /2 for any 9^. If true, this would mean that the asymptotic Higgs 
field configuration, $5, for the sphaleron for arbitrary 9^ is given by 

( sin^(0/2)e^^^+cos2(0/2) \ 
^ l^sin(6i/2)cos(6i/2) e''^(l-e-*''s)y' ^ ' 

where 75 = tt/ cos(20^). 

The twisting of the magnetic field lines in the sphaleron configuration has been further 
clarified in |Q . The direction of magnetic field lines is shown for a dumbell in Fig. |2^ and 
for a "stretched" sphaleron in Fig. ^ (The asymptotic fields for the stretched sphaleron 
are identical to those for the sphaleron and the twisted Z-string.) In the stretched sphaleron 
case, the magnetic field line twists around the vertical string segment by an angle tt (for 
Ow 0) as one goes from monopole to antimonopole. This twist provides non-trivial Chern- 



Simons number to the configuration |123| 



On physical grounds it seems reasonable that there should be a critical value of twist at 
which one can get a static solution for a Z-string segment. This is because the segment 
likes to shrink under its own tension but the twist prevents the shrinkage and is equivalent 
to a repulsive force between the monopole and antimonopole. (This idea owes its origin to 



Taubes |113| who discovered a solution containing a monopole and an antimonopole in an 
0(3) model in which the Coulomb attraction is balanced by the relative misorientation of 
the magnetic poles.) Then, if the string is sufficiently twisted, the attractive force due to the 
tension and the repulsive force due to the twist will balance and a static solution can exist. 
So far we have been assuming that the only dynamics of the segment is towards collapsing 
or expanding of the string segment. However, since we are dealing with twisted segments, 
we should also include the rotational dynamics associated with twisting and untwisting. So, 
while any twist greater than a certain critical twist might successfully prevent the segment 
from collapsing, only a special value of the twist can give a static solution to the rotational 
dynamics. Furthermore, we expect that this solution will be unstable towards rotations that 
twist and untwist the string segment. This would be the unstable mode of the sphaleron. 

Similar connections between the W-string and the sphaleron have also been constructed in 
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Figure 20: The thick sohd hne is the location of the Z-string for a dumbell configuration 
and the dashed curves he in the equatorial plane and are drawn to guide the eye. The 
dotted lines depict lines of magnetic flux. The arrows show the orientation of the vector 
n oc — i>^ri>. 
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Figure 21: The field configuration for a stretched sphaleron as in Fig. 
field line is shown. 



M Only one magnetic 




11 The ■^He analogy 



The symmetry structure of ^He closely resembles the electrowea k symmetry g roup and 



hence we expect the analog of electroweak strings to exist in "^He [132, 129, 131]. Indeed, 
this analog is called the n = 2 vortex. We now explain this correspondence in greater detail. 



11.1 Lightning review of ^ He 



•^He nuclei have spin 1/2 and two such nuclei form a Cooper pair which is the order parameter 
for the system. Unlike "^He, the pairing is a spin triplet (5 — 1) as well as an orbital angular 
momentum triplet (L = 1). As a result there are 3x3 components of the wavefunction of 
the Cooper pair - that is the order parameter has 9 complex components. Hence, the order 
parameter is written as a 3 by 3 complex valued matrix: Aai with a (spin index) and i 
(spatial index) ranging from 1 to 3. 

At temperatures higher than a few milli Kelvin the system is invariant under spatial ro- 
tations (50(3)^) as well as rotations of the spin degree of freedom of the Cooper pair 
(50(3)5). Another symmetry is under overall phase rotations of the wavefunction (U{1)n) 
and the corresponding conserved charge is particle number (N) . Hence the symmetry group 
is: 

G = SOi3)L X S0{3)s X U{1)n ■ (241) 
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There are several possible phases of ^He corresponding to different expectation values of 
the order parameter. In the A-phase, the orbital angular momenta of the Cooper pairs are 
all aligned and so are the spin directions. This corresponds to 

A^, = ^oL'4>^ (242) 

where Aq ~ 10~^ eV is the temperature dependent gap amplitude, the real unit vector da 
is the spin part of the order parameter, and 

tP, = — — — (243) 

with m and n being orthogonal unit vectors, is the orbital part of the order parameter. This 
expectation value of the order parameter leads to the symmetry breaking: 

G ^ [7(1)53 X U{1)l,-n/2 X Z2 . (244) 

The reason why a U{1) subgroup of SO{3)l x U{1)n survives the symmetry breaking can 
be derived from the expectation value in eq. ( |242| ) . A spatial rotation of the order parameter 
is equivalent to a phase rotation of ipi and this phase can be absorbed by a corresponding 
U{1)n rotation of the order parameter. Hence, just as in the electroweak case, a diagonal 
C/(l) subgroup remains unbroken. The U{l)s3 survives since rotations about the d axis leave 
the order parameter invariant. The non-trivial element of the residual discrete Z2 symmetry 
corresponds to a sign inversion of both tpi and d. A depiction of the A- and B- phases is 
shown in Fig. ^ (after [^). 

In the B-phase, neither the orbital angular momenta nor the spin directions of the different 
Cooper pairs are aligned. But the angle between the direction of the angular momenta 
and the spin direction is fixed throughout the sample. Hence in the B-phase, independent 
rotations of the orbital angular momenta and of spin are no longer symmetries. However, 
a simultaneous rotation of both orbital angular momenta and spin remains an unbroken 
symmetry. In other words, a diagonal subgroup of S0{3)s x SO{3)l remains unbroken. 
Therefore, in the B-phase the order parameter is written as: 

= 3-i/V'^i?„,(n,0) (245) 

where, is a phase and the 3x3 matrix Rai describes relative rotations of the spin and 
orbital degrees of freedom about an axis n and by angle 9. The symmetry breaking pattern 
is: 

G ^ SO{3)l,+Ss ■ (246) 
This symmetry breaking resembles the chiral symmetry breaking transition studied in QCD 
(with two flavors of quarks) and may be useful for experimentally investigating phenomenon 
such as the formation of "disoriented chiral condensates" |^ . The B-phase does not resem- 
ble the electroweak model and hence we will not discuss it any further. We shall also not 
discuss the various other phases of '^He (for example, the Ai phase) which are known to 
occur. (For a useful chart of the phases, see Sec. 6.2 of Ref. [ |132[ .) 

In addition to the continuous symmetries, there are a number of discrete symmetries that 
arise in the phases of ^He. These are important for the classification of topological defects 



in He. A description may be found in |111| 
11.2 Z-string analog in ^ He 



Clearly the A-phase closely resembles the electroweak symmetry breaking because of the 
mixing of the generator of the non-Abelian group (5*0(3)^) and the Abelian group (C/(l)jv)- 
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Figure 22: Depiction of the A- and B- phases of ^He. In the A-phase, the spin orientations 
of all the Cooper pairs are parallel and so are the orbital orientations. In the B-phase, the 
relative orientation of the spin and orbital orientations are fixed in all the Cooper pairs but 
neither the spin nor the orbital orientations of the various Cooper pairs are aligned. 




Spin orientation 



orbital orientation 



B-phase 
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The orbital part of the order parameter is responsible for this pattern of symmetry breaking 
and hence ipi plays the role of the electroweak Higgs field $. The connection, however, is 
indirect since tpi is a complex 3 vector while <J> is a complex doublet. The idea is that the 
■^He-A real vector 

IhbA = i-^—^-^ = rh X h (247) 



is analogous to the electroweak real vector 



$tf*(j) 



lero = -^TT-^ • (248) 

The electroweak Z-string is a non-topological solution for which the Higgs field configuration 
is: 

$ = -^/(r)e*'^(^5) . (249) 

For this configuration l^w = z. 



The vacuum manifold Ma of ''He-A has 



tti{Ma) = Z4 (250) 

and hence there are topological Z4 vortices in ■^He-A. The vortices occur in classes labeled 
by n = ±1/2, 1. The vortices with n equal to an even integer are topologically equivalent 
to the vacuum. The non-trivial topological vortices (labeled by n = ±1/2, 1) cannot be the 
equivalent of the non-topological Z-string. However, the topologically trivial n = —2 vortex 
is also seen in "^He-A. The order parameter for this vortex is: 

A^,{p, ^) = Aoz„[ e™'^/i(p)(i, + ly^) + e'^^'+^'^^ h{p){x, - ly,)] . (251) 

where fi{p) and f2{p) are two profile functions with /i(oo) = 1, /2(oo) — 0, /i(0) = and 
72(0) depending on n. In correspondence with the electroweak Z-string, the n = 2 vortex 
has IheA = z. However, the order parameter need not vanish at the center of the vortex for 
certain members of the n = 2 class of vortices. For example, with n = —2, we may have 
/2(0) ^ 0. 



The n — 2 vortex is not topological and can be continuously deformed into the vacuum 
manifold. The configuration at the terminus of the n — 2 vortex is called the hedgehog 
or monopole IhcA — f (the radial unit vector). This texture is the direct analog of the 
electroweak magnetic monopole {l^w = f) at the terminus of a Z-string. 

The n = 2 discontinuous vortex is unstable but even so has been observed in '^He. In the 
laboratory, the rotation of the sample stabilizes the n — 2 vortex. This seems to be closely 
analogous to the result of Garriga and MontesEq] who find that electroweak strings can 



be stabilized by external magnetic fields (Sec. 7.3) 



Before proceeding further, it is prudent to remind ourselves of some important differences 
between the (bosonic sector of the) GSW model and '^He. The symmetries in ^He are all 
global whereas the symmetries in the GSW model are all local. So the n = 2 discontinuous 
vortex is like a global analog of the Z-string. Another important difference is in the discrete 
symmetries in the two systems. The symmetry structure of the GSW model is really \SU (2) x 
U{l)]/Z2 since the Z2 elements 1 and —1 which form the center of SU (2) also occur in U (1). 
On the contrary, the symmetry group of '^He-A has a multiplicative Z2 factor which gives 
rise to the non-trivial topology of the vacuum manifold. 
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It is important to note that we cannot expect '^He to provide an exact replica of the GSW 
model. However, the similar structures of the two systems means that certain issues can be 
experimentally addressed in the '^He context while they are far beyond the reach of current 
particle physics experiments. An issue of this kind is the baryon number anomaly in the 
GSW model and the anomalous generation of momentum in ■^He. 

As described in Sec. ^, there are fermionic zero modes on the Z-string and an electric field 
applied along the Z-string leads to the anomalous production of baryon number. What is 
the corresponding analog in ■^He? At first sight, "^He does not have the non-Abelian gauge 
fields that the electroweak string has and so it seems that the analogy is doomed. But 
this is not true. The point is that the physics of fermionic zero modes has to do with 
the dynamics of fermions on the fixed background of the Z-string. Likewise, in "^He we 
can be interested in the dynamics of quasiparticles in the fixed background of the n = 2 
vortex. As far as the interaction of quasiparticles with the order parameter background 
is concerned, one can think of the '^He-A vortex as being due to a to a (fictitious) gauge 
field Z"*. Then the interaction of quasiparticles with the order parameter is of the form 
jf_iZ"^ which is exactly analogous to the interaction of quarks and leptons with the Z-boson. 
Just as in the electroweak case, the '^He quasiparticles have zero modes on the vortex. In 
close analogy with the scenario where the motion of a superconducting string through an 



external magnetic field leads to currents along the string (Sec. 3.4), the velocity of the "^He 
vortex through the superfluid leads to an anomalous flow of quasiparticles but this time 
in the direction perpendicular to the vortex. This flow causes an extra force on the vortex 
as it moves through the superfluid that can be monitored experimentally. Such a force 
was measured in the Manchester experiment JTo] , |13C| | and is in excellent agreement with 
theoretical predictions. Hence the Manchester experiment verifies the anomalous production 
of quasiparticle momentum on moving vortices and the corresponding production of baryon 
number on electroweak strings moving through a magnetic field. 



12 Concluding remarks and open problems 



Quantum field theory has been very successful in describing particle physics. Yet the suc- 
cesses have mostly been relegated to perturbative phenomena. A more spectacular level of 
success will be achieved when our field theoretic description of particle physics is confirmed 
at the non-perturbative level. The first non-perturbative objects that are likely to be en- 
countered in this quest are topological defects and their close cousins that we have described 
in this review. 

The search for topological defects can be conducted in accelerator experiments or in the 
cosmological realm via astronomical surveys. These searches are complementary - only su- 
permassive topological defects can be evident in astronomical surveys, while only the lightest 
defects can potentially be produced in accelerators. Foreseeable accelerator experiments give 
us access only to topological defects at the electroweak symmetry breaking scale. So it is 
very important to understand the defects present in the standard electroweak model and 
all its viable extensions. One may hope that the structure of defects will yield important 
clues about the underlying symmetry of the standard model. 

With this hope, we have described wide classes of defects present in field theories. These 
defects are not all topological and this is relevant to the standard electroweak model which 
also lacks the non-trivial topology needed to contain topological defects. The absence of 
topology in the model means that the defect solutions cannot be enumerated in topological 
terms and neither can their stability be guaranteed. We have described, however, how 
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the existence of defect solutions may still be derived by examining the topological defects 
occurring in subspaces of the model. The electroweak defects can be thought of as being 
topological defects that are embedded in the electroweak model. 

The issue of stability of the defect solution is yet more involved and has not yet been fully 
resolved in the presence of fermions. That the electroweak Z-string is stable for large 0^, 
(within the bosonic sector) was inspired by the discovery of semilocal strings and their 
stability properties. The explicit stability analysis of the electroweak string marks out the 
region of parameter space in which the Z-string is stable. Then it is clear that the Z-string 
is unstable for the parameters of the standard model. In certain viable extensions of the 
standard model and under some external conditions (such as an external magnetic field), 
the standard electroweak Z-string can still be stable. 

Even if the Z-string is unstable, it is possible that the lifetime of segments of string is long 
enough so that they can be observed in accelerators. This possibility was discussed in the 
first paper on the subject by Nambu p3| . The discovery of Z-string segments would truly 
be historic since it would confirm the existence of magnetic monopoles in particle physics. 
However, the rate of formation of Z-string segments and their lifetime has not yet been 
studied in detail. Some of the difficulties in this problem lie quite deep since they involve 
the connection of perturbative particle physics to the non-perturbative solitonic features. 
Additionally, the influence of fermions on electroweak strings needs further investigation. 

Electroweak strings may play a cosmological role in the genesis of matter over antimatter as 
is evident since configurations of electroweak string have properties that are similar to the 
electroweak sphaleron. The challenge here is to determine the number density of electroweak 
strings formed during the electroweak phase transition and their decay rate. Note that the 
formation of topological strings has been under constant examination over the last two 
decades and only now, with some experimental input, are we beginning to understand their 
formation. The cosmological formation of electroweak defects has not been addressed with 
as much vigour. Recently though, there have been spurts of activity in this area, with lattice 
calculations beginning to shed interesting insight p6[ . It is very likely that further lattice 
results will be able to give quantitative information about the formation of electroweak 
strings at the electroweak phase transition. 

While particle physics experiments to detect electroweak strings are quite distant, experi- 
ments in condensed matter systems to study topological defects are becoming more feasible 
and can be used to test theoretical ideas that are relevant to both particle physics and con- 
densed matter physics. Already there are experiments that test theories of the formation 
of topological vortices. We can also expect that condensed matter experiments might some 
day test the formation of defects that are not topological. The experiments on He"^ are most 
relevant in this regard since it contains close analogs of electroweak strings. Furthermore, 
ideas relating to the behaviour of fermions in the background of electroweak strings can 
also be tested in the realm of He"^. This makes for exciting physics in the years to come 
which will stimulate the growth of particle physics, cosmology and both, theoretical and 
experimental, condensed matter physics. 
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